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ABSTRACT 

We present Keck/LRIS spectra of over 200 galaxies with well-determined redshifts between 0.4 and 
1.4. We combine new measurements of near-ultraviolet, low-ionization absorption lines with previously 
measured masses, luminosities, colors, and star formation rates to describe the demographics and 
properties of galactic flows. Among star-forming galaxies with blue colors, we find a net blueshift 
of the Fe II absorption greater than 200 km s _1 (100 km s _1 ) towards 2.5% (20%) of the galaxies. 
The fraction of blueshifted spectra does not vary significantly with stellar mass, color, or luminosity 
but does decline at specific star formation rates less than roughly 0.8 Gyr -1 . The insensitivity of 
the blueshifted fraction to galaxy properties requires collimated outflows at these redshifts, while 
the decline in outflow fraction with increasing blueshift might reflect the angular dependence of the 
outflow velocity. The low detection rate of infalling gas, 3 to 6% of the spectra, suggests an origin in 
(enriched) streams favorably aligned with our sightline. We find 4 of these 9 infalling streams have 
projected velocities commensurate with the kinematics of an extended disk or satellite galaxy. The 
strength of the Mg II absorption increases with stellar mass, B-band luminosity, and U-B color, trends 
arising from a combination of more interstellar absorption at the systemic velocity and less emission 
filling in more massive galaxies. Our results provides a new quantitative understanding of gas flows 
between galaxies and the circumgalactic medium over a critical period in galaxy evolution. 

Subject headings: 



1. INTRODUCTION 

Key processes in the evolution of galaxies include the 
flow of cool gas into galaxies, the conversion of these 
baryons into stars, and the ejection of gas enriched with 
heavy elements. Determining the factors that govern the 
circulation of baryons remains a critical, unsolved prob- 
lem in cosmology. 

The measured properties of galaxies indicate strong 
evolution in this baryon cycle between the peak in star- 
formation activity at redshift z k, 2 and the current 
epoch (Lilly etal. 1996; Madau etal. 1996; Hopkins & 
Beacom 2006). At any redshift during this period, star- 
forming galaxies populate a fairly well-defined locus in 
the star formation rate (SFR) - stellar mass plane (Bell 
etal. 2005; Elbaz etal. 2007; Noeske etal. 2007). How- 
ever, the observed SFR at a fixed stellar mass declines by 
about a factor of 20 from z ~ 2 to the present. The Tully- 
Fisher relation also evolves in the sense that present-day 
galaxies have about 2.5 times more stellar mass than do 
galaxies at z 2.2 with the same circular velocity (Cresci 
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etal. 2009). Since redshift z ~ 1 when roughly half the 
mass in present-day galaxies was assembled (Drory et al. 
2005; Faber et al. 2007; Marchesini et al. 2009), the stellar 
mass in red sequence galaxies has continued to grow; but 
the mass in blue galaxies has remained essentially con- 
stant (Faber etal. 2007). Over this period, some process 
quenches star formation in massive galaxies, yet observa- 
tional efforts to identify this process with feedback from 
active galactic nuclei (AGN) have not reached a consen- 
sus (cf., Schawinski etal. 2007; Aird etal. 2012). 

The leading explanation for this drop in cosmic star- 
formation activity is a decline in the accretion rate of 
cool gas onto galaxies. The simplest form of this idea, 
that the mean growth rate of dark matter halos regulates 
the gas inflow rate onto galaxies, fails; it yields too much 
star formation in low mass halos at early times (Bouche 
et al. 2010) and does not shut down star formation in high 
mass halos at later times (Somerville etal. 2008). Intro- 
ducing gas physics - heating by virial shocks in massive 
halos and outflows in both low-mass and high-mass ha- 
los - effectively tilts the underlying relationship between 
halo accretion rate and mass (described by the simplest 
model) into the observed SFR - stellar mass relation and 
Tully-Fisher relation (Bouche et al. 2010 ; Dave, Oppen- 
heimer, & Finlator et al. 2011a; van de Voort et al. 2011a, 
b; Dave, Finlator, & Oppenheimer 2012). The balance 
between gas outflow and inflow then determines the re- 
lationship between the gas-phase metallicity and stel- 
lar mass (Dave, Finlator, &: Oppenheimer 2011b). The 
problem with this scenario is that this baryon-driven pic- 
ture of galaxy formation rests on a rather rudimentary 
understanding of which galaxies have outflows and in- 
flows and how the physical properties of flows vary with 
fundamental galaxy properties. 
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Gas flows imprint resonance absorption lines on the 
spectra of their host galaxies that can often be kinemat- 
ically distinguished from interstellar gas at the systemic 
velocity (Heckman et al. 2000; Schwartz & Martin 2004; 
Martin 2005, 2006; Rupke etal. 2005; Schwartz ctal. 
2006; Tremonti et al. 2007; Martin & Bouche 2009; Rubin 
et al. 2010a). Thus far, studies of gas flows at 0.5 < z < 2 
have either been based on just a few individual galaxies 
(Sato etal. 2009; Coil etal. 2011; Rubin etal. 2011) or 
on composite spectra, i.e., the average of many low S/N 
ratio spectra (Weiner etal. 2009; Rubin etal. 2010b). 
The pioneering study of Weiner et al. (2009) concluded 
that blueshifted Mg II absorption was ubiquitous in spec- 
tra of star- forming galaxies at z <~ 1.4 and demonstrated 
that the Doppler shift and absorption equivalent width of 
this absorption increase (rather slowly) with both stellar 
mass and SFR. At lower redshift, however, Rubin et al. 
(2010) did not find blueshifted Mg II absorption in a com- 
posite spectrum of galaxies with stellar masses similar 
to the Weiner etal. (2009) sample and speculated that 
the higher specific SFRs (i.e., SFR per unit stellar mass) 
characterizing z « 1.4 galaxies might be required to host 
such outflows. These two studies have demonstrated the 
existence of outflows and their possible evolution over a 
key period in the assembly of galaxies but are limited 
by their reliance on the mean spectrum of a population. 
Averaging many spectra together in a composite hides 
less common features like gas inflows (Sato et al. 2009; 
Rubin etal. 2012) and does not allow de-projection of 
the line-of-sight Doppler shift into an outflow velocity. 

Here, we present the results of a survey of near- 
ultraviolet spectral features in 208 galaxies with redshifts 
between 0.4 < z < 1.4 to provide an empirical mea- 
surement of how outflow and inflow properties change 
with cosmic time and galaxy properties. The galaxies 
were selected from the DEEP2 (Deep Extragalactic Evo- 
lutionary Probe 2) survey (Davis etal. 2003; Newman 
et al. 2012). This redshift survey provides a relatively un- 
biased sample for investigating the demographics of out- 
flows and inflows in galaxies brighter than Rab = 24.1 
in four fields. More importantly though, selecting from 
the DEEP2 survey ensures that fundamental galaxy pa- 
rameters - such as stellar mass, i?-band luminosity, and 
U — B color - have been measured for the sample in a sys- 
tematic manner (Bundy etal. 2006; Willmer etal. 2006). 
The additional photometry obtained for the AEGIS (All- 
Wavelength Extended Groth Strip International Survey) 
field allows us to measure SFRs for 51 of these galaxies; 
and these SFRs provide important information about the 
energy and momentum produced by supernovae, stellar 
winds, and radiation from massive stars and therefore 
available to drive outflows (Chevalier & Clegg 1985; Mur- 
ray et al. 2005). 

In this paper, we focus on the diagnostics provided 
by low-ionization, resonance-absorption lines, reserving 
the presentation of high-ionization interstellar absorp- 
tion and stellar features, resonance emission, and fluo- 
rescent emission for future papers. Low-ionization, res- 
onance lines in near-ultraviolet (rest-frame) spectra pro- 
vide an empirical bridge between the optical transitions 
typically observed for low-redshift galaxies and the far- 
UV transitions studied extensively in spectra of high- 
rcdshift galaxies (Steidel et al. 2010) because the Mg II 
AA2796, 2803 doublet is accessible from the ground over 



the broad redshift range from roughly 0.25 < z < 2.5 
(Martin & Bouche 2009; Tremonti etal. 2007; Weiner 
etal. 2009; Rubin etal. 2010; Coil etal. 2011). At the 
lower end of this redshift range, Mg II absorption prop- 
erties can be cross calibrated with rest-frame optical lines 
such as Na I AA5890, 96 and Ca II A3933, 69 . At z> 1.19, 
the Mg II absorption properties can be directly compared 
to far-UV transitions in Si II, Al II, C II, and C IV 
commonly employed to study outflows in much higher 
redshift galaxies (Shapley etal. 2003; Steidel etal. 2010; 
Jones etal. 2012). 

The Mg II doublet provides a very sensitive probe of 
outflows for several reasons: singly-ionized magnesium is 
a dominant ionization state over a broad range of condi- 
tions, both lines have large oscillator strength, and the 
cosmic abundance of Mg is fairly high. Scattered Mg II 
emission, however, partially fills in the intrinsic absorp- 
tion troughs in galaxy spectra and complicates the inter- 
pretation of the Doppler shift of the absorption trough. 
Spectral coverage blueward of the Mg II AA2796, 2803 
doublet provides access to a series of strong Fe II reso- 
nance lines which alleviate the concern about emission 
filling. The Fe II absorption troughs provide a cleaner 
view of the intrinsic absorption profile because fluores- 
cence (rather than resonance emission) often follows ab- 
sorption in several of these Fe II transitions (Prochaska 
etal. 2011; Erb etal. 2012). Another advantage of Fe II 
lines over the Mg II doublet arises from the large number 
of Fe II transitions in the NUV. The oscillator strengths 
of the NUV Fe II transitions span a substantial range and 
therefore make it possible to place useful bounds on the 
column density of singly-ionized iron, thereby constrain- 
ing the total gas columns for an assumed (i.e., model 
dependent) metallicity and ionization fraction. 

The paper is organized as follows. Section 2 introduces 
the sample and describes our new Keck observations, in- 
cluding a discussion of the systemic velocity determina- 
tion and absorption-line sensitivity. After providing a 
broad overview of the ncar-UV spectral features, Sec- 
tion 3 explains the complications caused by emission fill- 
ing and how we identify and measure galactic gas flows. 
We then compare the line profiles of the Mg II and Fe II 
absorption troughs in Section 4. In Section 5, we cal- 
culate the fraction of galaxies with a Doppler shift of 
low-ionization absorption relative to the systemic ve- 
locity and use previously measured galaxy properties 
to illustrate the demographics of the galaxies showing 
blueshifts. We then discuss the physical properties of 
the outflows and show how the outflow properties scale 
with galaxy properties. In Section 6, we expound on 
our discovery of net inflows of enriched gas. Section 7 
summarizes our conclusions. 

Throughout this paper we assume a cosmological 
model with tt m — 0.3, fl\ = 0.7, and Hq 
71 1 km s _1 Mpc -1 . We adopt the atomic oscillator 
strengths and cosmic abundance ratios given by Mor- 
ton (2003) as well as the associated vacuum wavelengths 
for transitions shortward of 3200A. We refer to optical 
transitions by their wavelengths in air (for ease of com- 
parison to previously published work) , but we work with 
their vacuum wavelengths. 
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2. DATA 

We present rest-frame, ultraviolet spectroscopy of z <~ 

I galaxies drawn from the Deep Extragalactic Evolution- 
ary Probe 2 survey (Davis etal. 2003; Newman etal. 
2012) using the Low Resolution Imager and Spectrometer 
(LRIS) on the Keck I telescope (Oke etal. 1995). Build- 
ing a sample of 208 spectra required a significant obser- 
vational campaign, custom-designed masks for the focal 
plane, and the remarkable blue sensitivity of Keck/LRIS- 
B (Steidel etal. 2004). We describe the sample selection 
in Section 2.1 and the new observations in Section 2.2. 

2.1. Sample 

The DEEP2 survey obtained redshifts for approxi- 
mately 38,000 galaxies spread over four fields widely 
separated in right ascension (Newman etal. 2012). In 
three of the four fields, the galaxies were color-selected 
to have redshifts z > 0.7. The fourth field, without 
color selection, is also the target of a deep imaging ef- 
fort from the X-ray through the radio known as the All- 
Wavelength Extended Groth Strip International Survey, 
AEGIS (Davis etal. 2007). The panchromatic data in 
AEGIS provide more accurate estimates of galactic SFRs 
than are currently available for the other three fields 
(e.g., Noeske etal. 2007). Figure 1 summarizes the prop- 
erties of the galaxies observed with LRIS. 

Our primary sample of galaxies was chosen for spectro- 
scopic follow-up based on redshift, 1.19 < z < 1.35, and 
apparent magnitude B < 24.0. The magnitude cut en- 
sures adequate continuum S/N ratio for absorption-line 
spectroscopy. This lower bound on the redshift guar- 
antees that the C IV AA1548, 1550 doublet lands long- 
ward of the atmospheric cut-off in the blue; and requiring 
z < 1.35 ensures that the Mg I 2853 transition is covered 
by the blue spectrum. The simultaneous coverage of far- 
ultraviolet (FUV) and near-ultraviolet (NUV) resonance 
lines allows a direct comparison of the transitions used 
to measure outflow velocities in z > 2 galaxies (Steidel 
etal. 2010) and z <~ 1 galaxies (Weiner etal. 2009; Rubin 
etal. 2010). In total, 68 galaxies meeting these criteria 
were observed on 6 masks. Each mask subtends a 5.'5 by 
8'.0 field of view (Steidel etal. 2004). With an average of 

II primary targets per mask, we filled the masks with the 
following types of galaxies: (1) lower redshift (z < 1.19) 
galaxies with EW([ O II ]) > 10 A and B < 24.5 mag- 
nitudes, green valley galaxies, and K+A galaxy candi- 
dates (Yan etal. 2009). In total, 145 galaxy spectra were 
obtained through these masks with the 400 lines mm -1 
configuration of LRIS described below. 

To utilize less than optimal observing conditions, a sec- 
ond set of masks were designed to obtain NUV spec- 
troscopy of brighter galaxies. Spectra were obtained 
through 3 masks with galaxies prioritized as follows: (1) 
B < 23.6, z > 0.7, and EW([ O II ]) > 20 A (41 ob- 
jects); (2) B < 23.6, z > 0.7, and EW{[ O II ]) < 20 A 
(12 objects); (3) B = 23.6 - 24.0 and z > 0.7 (7 objects), 
B < 23.6 and z < 0.7 ( 1 object), and 6 foreground 
galaxies that lie at small angular separation from an- 
other target. Among these 67 galaxies, the spectra of 
four foreground objects with z < 0.4 are not discussed 
further in this paper due to their much lower redshifts. 
Hence, the second sample includes 63 galaxies observed 
with the 600 lines mm -1 configuration of LRIS described 



below. 

The primary sample (z > 1.19) includes three galaxies 
with colors that place them between the red sequence 
and blue cloud, a region of the color-magnitude diagram 
known as the green valley (Bell et al. 2003; Mendez et al. 
2011). A total of 21 green valley galaxies were included 
on the secondary (600 lines mm -1 backup) and primary 
(400 lines mm" 1 ) masks including the filler objects; most 
have EW([OII]) < 10 A. The sample includes just two 
K+A (post-starburst) galaxies identified by Yan et al. 
(2009). One of these is a green valley galaxy and the 
other is a K+A galaxy with colors just redward of the 
green valley in the color - magnitude diagram. 

As illustrated in Figure 1, these targets span the red- 
shift range from 0.4 to 1.4, have luminosities between 
— 17 > Mb > —22, and sample the stellar mass func- 
tion between 8.9 < log(M/ M Q ) < 11.6. With a few 
exceptions noted above, the galaxies have blue U — B 
colors. These spectra therefore probe the properties of 
low-ionization gas in and around typical star-forming 
galaxies 4.3 to 9.0 Gyr ago. The broad range in stel- 
lar mass makes this an appropriate sample to compare 
gas outflows and inflows over about 3 dex in dark-matter 
halo mass. The clustering of the population studied sug- 
gests that the most massive fertile of our sample re- 
sides in halos with log(M/j/ M Q ) > 12 and that some 
of the galaxies in the lowest mass fertile populate ha- 
los with log(M/j/ M Q ) < 11. In the stellar mass range 
log(M*/ M Q ) = 9.0 — 11.5, abundance matching sug- 
gests \og(M h / M Q ) = 11.25 - 14.24 at z = 1.0 and 
log(AV M ) = 11.2 - 14.7 at z = 0.5 (Behroozi, Con- 
roy, and Wechsler 2010), consistent with the halo masses 
measured from galaxy clustering. 

2.2. LRIS Observations 

Multislit spectroscopy of the fields listed in Table 1 was 
obtained with LRIS on the Keck I telescope between 2007 
October and 2009 June. Our primary instrumental con- 
figuration used the D680 dichroic, the LRISb 400 mm" 1 
grism blazed at A3400A, and the the 831 mm -1 grating 
blazed at 8200A on the red side. For the backup program, 
we configured LRIS with the D560 dichroic, 600 mm -1 
grism blazed at 4000A, and the 600 mm red grating 
blazed at 7500A; this setup generally provided continu- 
ous spectral coverage between the red and blue spectra. 
In both configurations, the red-channel wavelength cov- 
erage usually included rest-frame-optical emission lines, 
most commonly [O II] AA3726, 29. For all observations, 
the Cassegrain Atmospheric Dispersion Corrector (ADC) 
(Phillips et al. 2006) was enabled; and the width of the 
slitlets subtended l'/2 on the sky. We observed a total of 
9 slitmasks, 6 in our primary configuration, and 3 in our 
back-up configuration, as described below. Conditions 
during the observing runs ranged from excellent (clear 
with 0'.'6 seeing) to variable (intermittent clouds with 
variable seeing up to 2'.'0). Details of these observations 
are presented in Table 1. 

The data were primarily reduced using IRAF tasks, 
with scripts designed for cutting up the multi-object slit 
mask images into individual slitlets, flat-fielding using 
spectra of the twilight sky for the blue side and dome 
flats for the red side, rejecting cosmic rays, subtracting 
the sky background, averaging individual exposures into 
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FlG. 1. — Properties of the galaxies observed with LRIS. (Top:) Redshift histogram. Higher redshift galaxies were specifically targeted 
for the lower resolution observations in order to obtain spectral coverage of resonance lines in the far-UV spectrum. The median redshifts of 
the galaxies observed with the 600 lines mm -1 (282 km s _1 FWHM) and 400 lines mm -1 (435 km s _1 FWHM) configurations are 0.94 
and 1.12, respectively. (Middle:) Color - magnitude diagram (Willmer etal. 2006). The diagonal, dashed line marks the division between 
the red sequence and the blue cloud at z ~ 1. For the LRIS sample, galaxy color is anti-correlated with absolute magnitude, 7.8 standard 
deviations from the null hypothesis (i.e., no correlation) with Spearman rank correlation coefficient r$ = —0.54. The rest-frame optical 
colors of the galaxies (Willmer et al. 2006) become slightly bluer with increasing redshift due to the .R-band selection of DEEP2 galaxies, and 
the lowest redshift fertile includes lower luminosity galaxies that would not pass our apparent magnitude cut at higher redshift. (Bottom:) 
Color vs. stellar mass (Bundy et al. 2006) showing histogram of stellar mass along the bottom. Stellar mass was derived from SED fitting 
using a Chabrier stellar initial mass function (IMF; Chabrier 2003). The stellar masses would be 0.25 dex larger for a Salpeter IMF. The 
histogram distinguishes the most secure stellar masses (solid bars), obtained for 127 galaxies detected in the K-band, from the full sample 
of 208 galaxies with LRIS spectroscopy. The median stellar mass of the K-band-detected galaxies in our sample is 1.60 X 10 10 Mq, and 
the galaxies not detected in the K-band have a median mass of 7.23 X 10 9 Mq as estimated from optical photometry. Color and stellar 
mass are strongly correlated (rg = 0.77 at llcr significance). 



final stacked two-dimensional spectra, extracting to one 
dimension, wavelength and flux calibrating, and shifting 
into the vacuum frame. 10 These procedures are described 
in detail in Steidel et al. (2003). We also followed the pro- 
cedures outlined in Shapley et al. ( 2006) for background 

10 The typical root-mean-square error in the wavelength calibra- 
tions are 0.30 A, 0.15 A, 0.05 A, and 0.05 A, respectively, for the 
blue 400 lines mm -1 , blue 600 lines mm -1 , red 831 mm -1 , and 
red 600 mm -1 spectra. 



subtraction of deep mask spectra. Accordingly, to avoid 
potential oversubtraction of the background, the object 
continuum location was excluded from the estimate of 
the background fit at each dispersion point. Along with 
each target science spectrum, we extracted an error spec- 
trum from the standard deviation of the mean of the 
individual exposures. We paid particularly close atten- 
tion to deriving accurate wavelength solutions, adjust- 
ing an initial wavelength calibration derived from arc 
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lamp exposures such that sky lines in our spectra ap- 
peared at the correct wavelengths. Finally, spectra were 
flux calibrated using observations of spectrophotometric 
standards throughout the run. To calibrate the entire 
observed wavelength range, spectra of these stars were 
obtained through slits at multiple positions in the focal 
plane because the standard longslit spectral coverage is 
too blue relative to that of the multislit coverage. 

The spectral resolution realized differed between the 
two configurations. In the primary configuration 
(400 lines mm -1 grism), the atmospheric seeing rather 
than the slit width determined the spectral resolution 
because the angular sizes of the galaxies (measured along 
the spatial direction of the slit) were slightly smaller than 
the slit. Scaling the width of arc lamp lines by this ratio, 
the average effective resolution was 435 km s _1 FWHM 
(full width at half maximum intensity) in the grism 
spectra and 150 km s _1 in the grating spectra. In 
the 600 lines mm -1 configuration (used to observe lower 
redshift galaxies in non-optimal conditions), the galax- 
ies typically filled the l'/2 aperture; and the slit width 
set the spectral resolution. We measured line widths 
of 282 km s" 1 FWHM in the blue spectra and 220 
km s _1 FWHM in the red spectra. 

2.2.1. Determination of the Systemic Velocity 

Identifying net outflows or inflows requires accu- 
rate relative redshifts between the galaxy and the low- 
ionization gas. We derive the galaxy redshifts directly 
from the LRIS spectra whenever possible. The ADC 
used for the LRIS observations compensates for atmo- 
spheric differential refraction and eliminates systematic 
errors caused by the aperture subtending different re- 
gions of a spatially-resolved galaxy at widely separated 
wavelengths. In most of our spectra, line emission in 
the [O II] doublet drives the determination of the galaxy 
redshift. Since dense, photoionized gas near star-forming 
regions produce this emission, we expect these redshifts 
to accurately describe the systemic velocity. 

We derived LRIS redshifts for 167 galaxies using only 
the spectrum redward of the strong NUV Mg II and Mg I 
transitions (A r > 3000A). The most prominent spectral 
feature is usually the [O II] AA3726,29 doublet. How- 
ever, strong [O III] AA4959, 5007 and Balmer emission 
lines appear in spectra of the lower redshift galaxies; and 
Balmer absorption and the 4000A break are prominent 
in some spectra. Although we have attempted to avoid 
strong resonance lines that are blends of photosphcric 
and interstellar lines, interstellar Ca II A3933, 69 absorp- 
tion could lower some of the LRIS redshifts relative to the 
purely stellar/nebular redshift in a few spectra. To esti- 
mate a redshift, each spectrum was cross-correlated with 
3 templates - an emission line galaxy, a quiescent galaxy, 
and a K+A galaxy - using the DEEP2 IDL pipeline; the 
best fit to a linear combination of these 3 templates was 
adopted as the LRIS redshift. How accurately we can 
find the centroid of spectral lines determines the statisti- 
cal error on the redshift (and systemic velocity) for most 
galaxies. For a Gaussian line profile, the standard devi- 
ation of the mean, SDOM w a/SNR describes the cen- 
troiding error. For a single line, the typical uncertainty 
is 5{v) w 19 km s" 1 (cr/190 km s" 1 )(5/iV/10)- 1 , where 
a = FWHM (km s _1 )/2.35 is the spectral resolution. 



For the 41 galaxies with no LRIS coverage of [O II], we 
adopted the DEEP2 redshift and estimated the typical 
error in the systemic velocity from the redshift differ- 
ences measured among the 167 galaxies. A positive red- 
shift difference, zueepi — zlris > 0, means that adopt- 
ing the DEEP2 redshift increases the inferred outflow 
speed (makes it more negative) by an amount Av sys — 
c(z DE EP2 - zlris) /(I + zlris)- The mean offset com- 
puted from 167 spectra is Av sys = —14 ± 3 km s _1 , 
quite small compared to typical outflow/inflow speeds. 
The standard deviation of 41 km s _1 indicates the size 
of the systematic error in the redshift estimate. For a 
spatially-resolved source, uncertainties of this magnitude 
can result from an offset between the slit center and the 
brightest line emission; and although we use the same 
galaxy coordinates as the DEEP2 survey, the position 
angle of the slitlcts generally differ. 

2.2.2. Sensitivity to Absorption Lines 

Prior to making measurements, we deredshifted the 
spectra to the rest-frame. We fit a low-order spline 
to the continuum in each galaxy spectrum avoiding re- 
gions near spectral lines. We divided the spectra by the 
fitted continuum to produce normalized spectra. Un- 
certainty about the continuum level propagates into an 
error term in the equivalent width measurements. We 
found that the error term from the continuum fit could 
be kept small, compared the term derived from the error 
spectrum, by fitting the NUV continuum (approximately 
A2200-2900) independently of the FUV and optical con- 
tinua. 

The resulting sensitivity to absorption lines is de- 
scribed by the minimum detectable equivalent width. 
Our observing strategy was designed to reach a rest- 
frame equivalent width W r (2796) ~ 1 A because the 
distribution function of intervening absorbers presents 
a break here which may signal association with galaxy 
halos (Nestor etal. 2006). The weakest line that we 
detect depends on the continuum S/N ratio at the ob- 
served wavelength of the transition as well as the intrinsic 
width of the absorption trough. In the 2400-2500A band- 
pass, the median S/N ratios of the 400 lines mm -1 and 
600 lines mm - spectra are 6.5 and 5.5 (per pixel), re- 
spectively. For purposes of illustration, a typical equiv- 
alent width limit can be estimated assuming S/N ~ 5 
and the FWHM of an unresolved line. In the near-UV 
Fe II series, absorption troughs with rest-frame equiva- 
lent widths stronger than 

W 3a = 0.92 A {6.5/SNR) X2450 (1) 

in 400 lines mm -1 spectra, or larger than 

W 3a = 0.65 A {5.0/SNR) X2450 (2) 

in 600 lines mm -1 spectra, are easily detected. Just 
blueward of the Mg II doublet, the median continuum 
S/N ratios of the 400 lines mm -1 and 600 lines mm -1 
spectra are, respectively, 4.6 and 6.6. The typi- 
cal rest-frame sensitivities of the 400 lines mm -1 and 
600 lines mm - spectra are, respectively, 

W 3a = 2.02 A (5.0/SAi?) A2790 (3) 

and 

W 3a = 0.58 A (5.0/SNR) X2790 . (4) 
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Farther to the red near the Mg II doublet, the contin- 
uum S/N in a particular spectrum is usually a bit lower 
than it is near Fe II AA2587, 2600 due to the blue blaze 
of both grisms, so some of the spectra do not quite reach 
the target sensitivity in Mg II. Spectra with S/N ratio 
adequate to detect equivalent withs of 1-2 A in individual 
galaxies at z <~ 1 have never been previously presented 
however in such large numbers. 

We measure the Doppler shift of resonance absorp- 
tion relative to the systemic velocity. Since this sys- 
temic velocity, as well any projected gas flow, may vary 
slightly with slit position, the systematic error in the 
Doppler shift of the resonance lines could, in excep- 
tional geometries, be as large as the systematic error of 
41 km s" 1 estimated above. Typically, however, the 
statistical error in centroiding the absorption lines will 
characterize the error on the Doppler shift. The reso- 
lution of the 400 lines mm -1 spectra, expressed by the 
standard deviation a, is 1.54 times coarser than that of 
the 600 lines mrn^ 1 spectra; but this lower resolution is 
partly offset by higher average S/N ratio spectra of the 
400 lines mm -1 subsample. For emission lines in our 
best 400 lines mm -1 and 600 lines mm -1 spectra, this 
statistical error is about SV w 185/25 w 7 km s" 1 and 
5V « 120/10 s=s 10 km s _1 , respectively. Fitting mul- 
tiple absorption lines, as described in Section 3.3, fur- 
ther reduces the statistical error on our best estimate of 
Doppler shift V\ of the Fe II series. 

3. DIAGNOSTICS OF OUTFLOWS AND INFLOWS 

In this paper, we explore the diagnostics in the near- 
UV absorption-line spectrum longward of 2000 A. The 
emission line spectra are discussed further in two com- 
panion papers (Martin et al. , in prep; Kornei et al. , in 
prep). We defer the discussion of the far-UV spectral 
features to a future paper. 

3.1. Near-UV Spectral Features 

In Figure 2, we show the spectral region 
around the following NUV resonance lines: Fe II 
A2249.88, 2260.78, 2344.21, 2374.46, 2382.77, 2586.65, 
and 2600.17; Mg II AA2796.35, 2803.53; and Mg I 
A2852.96. The absence of strong absorption in some 
individual spectra will partially fill in the absorption 
troughs in the composite spectra. Inspection of the indi- 
vidual spectra, however, shows that other factors prevent 
the absorption troughs from appearing completely black. 
Foremost, instrumental resolution smooths the spectrum 
and explains why the absorption troughs are slightly 
deeper in the 600 lines mm -1 composite than in the 
400 lines mm -1 composite. In addition, the absorbing 
gas likely only partially covers the continuum source 
allowing some continuum light to leak out and fill in 
the absorption troughs; this partial covering has been 
fitted as a function of velocity using higher resolution 
spectra of brighter galaxies (Martin & Bouche 2009). 
All three of these factors - partial covering, instrumental 
resolution, and averaging spectra - imply that saturated 
absorption troughs need not be black. 

Nebular emission lines mark dense gas near star- 
forming regions excited by either photoionization or 
shocks. The LRIS red-side spectra typically show 
strong [O II] AA3726, 29 emission and cover an increas- 



ing number of rest-frame optical lines with decreasing 
redshift. In the NUV composite spectra, the [O II] 
AA2470.97, 2471.09 lines are detected and blended. We 
attribute the broader width of the strong emission fea- 
ture at A2326 to a blend of various C II] transitions; in 
our picture, the Fe II A2328 feature, which Rubin et al. 
(2010a) suggested contributed to this emission complex, 
does not contribute because its upper level cannot be 
populated by a permitted transition from the ground 
state. 

Line blending complicates measurements of the absorp- 
tion troughs. The individual transitions of the Mg II dou- 
blet blend together in the 600 lines mm^ 1 composite in 
Figure 2. The Fe II AA2587, 2600 doublet can blend with 
the Mn II AAA2576.88, 2594.50, and 2606.26 triplet, rec- 
ognizable in the higher resolution composite in Figure 2. 
Strong CII] A2326 emission often blends with the blue 
wing of A2344. In some individual spectra, the Mn II 
or CII] contamination compromises measurements of the 
maximum blueshift of absorption in the Fe II A2587 and 
Fe II A2344 line profiles, respectively. In contrast, the 
Fe IFA2365 emission blueward of Fe II A2374 is much 
weaker than the CII] A2326 emission, leaving A2374 as 
the preferred transition for looking for highly blucshiftcd 
Fe II absorption. 

The Fe II A2374 transition has the lowest oscillator 
strength among the NUV Fe II transitions routinely de- 
tected in the individual spectra. The oscillator strength 
of Fe II A2374 is only a tenth that of the strongest transi- 
tion, Fe II A2383. The oscillator strengths of Fe II A2383 
and Mg II A2803, the weaker line in the Mg II doublet, 
are nearly equal. The shorter wavelength transition in 
the Mg II doublet has an oscillator strength twice that 
of the longer wavelength transition. The rough equality 
of the absorption equivalent widths in all these transi- 
tions in Figure 2 demonstrates the resonance absorption 
is typically optically thick. The weak Fe II A2250 and 
A2261 transitions, just barely detected in Figure 2, are 
optically thin and particularly useful for constraining the 
column density of single-ionized iron, N(Fe + ). 

Among our spectra, we often find suprisingly low Fe II 
A2383 equivalent widths relative to the other Fe II ab- 
sorption lines. We attribute this apparent contradic- 
tion of the intrinsic line strengths to resonance emission 
partially filling in the deeper, intrinsic absorption pro- 
file. This emission filling is more significant in Fe II 
A2383 than the other strong Fe II transitions because 
Fe II A2383 absorption does not produce fluorescent emis- 
sion; i.e., the only permitted transition from the excited 
state is decay to the ground state. When an absorp- 
tion trough is partially filled in by emission, the centroid 
will generally differ from the instrinic absorption pro- 
file because absorption samples only the gas between the 
observer and the galaxy, while emission may come from 
both the near and far sides of the galaxy. 

Many of the individual spectra show prominent Mg II 
emission lines. In Figure 2, the longer wavelength, lower 
oscillator strength transition in the Mg II doublet has 
the larger equivalent width, an unphysical situation for 
the intrinsic equivalent width ratio. The inverted ap- 
parent ratio may arise from the troughs being partially 
filled in by resonance emission. For example, if the in- 
trinsic absorption troughs are saturated, then the two 
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Fig. 2. — Composite, near-ultraviolct spectra of star-forming galaxies at z ~ 1. The continuum level has been normalized to unity. The 
lower spectrum in Figure 2 shows the addition of all 208 continuum normalized spectra after smoothing them to a common resolution of 
435 km s _1 FWHM. The upper spectrum is the average of 63 spectra observed at a resolution of 282 km s —1 FWHM. Vertical blue (solid), 
red (dashed), and green (dotted) lines mark, respectively, resonance absorption lines, fluorescent emission lines, and nebular emission lines. 
The unlabeled Fe II* transition is a weak line at A2381.49; see Martin et al. (2012b, in prep.) for a tabular linelist. 



troughs will have equal intrinsic equivalent width; and 
slightly stronger emission in Mg II 2796 will produce 
the inverted ratio of the net equivalent widths. Simple, 
radiative-transfer models do not produce these inverted 
ratios (Prochaska etal. 2011), and more detailed model- 
ing is needed to test this idea. A non-spherical geometry, 
for example, allows the absorption and emission to come 
from physically distinct regions of the galaxy with differ- 
ing optical depths. Furthermore, it may be quite impor- 
tant to consider multiple origins of the Mg II line pho- 
tons, which are both emitted by HII regions (Erb et al. 
2012) and generated by continuum absorption through- 
out the interstellar medium (ISM) and circumgalactic 
medium. 

3.2. Emission Filling 

In the normal situation for intervening absorption in 
quasar spectra, where the continuum source subtends a 
small angular size compared to the gas clouds, a negligi- 
ble fraction of photons absorbed in the resonance tran- 
sitions are re-emitted along our sightline. For extended 
sources like galaxies, however, the equivalent width of the 
emission can be comparable to the absorption equivalent 
width. 

At low densities, every excitation from the ground state 
decays as either resonance emission or fluorescent emis- 
sion. In the absence of dust, photon number is conserved, 
and the equivalent width of the resonance emission (plus 
any fluorescent emission from the same upper level) will 
be exactly equal to the absorption equivalent width. The 
resonance emission will generally not cancel out the ab- 
sorption trough due to their velocity offset, generated 
for example by radial gas flow. The absorption trough 
reflects only the kinematics of gas in front of the con- 
tinuum source while the emission-line profile probes the 
kinematics of both the near and far sides of the galaxy. 



Departure from a spherically symmetric gas distribution 
generates specific viewing angles for which spectra will 
show stronger absorption than emission (and other an- 
gles with emission stronger than absorption). In addition 
to dust and the geometry of the gas distribution, spec- 
troscopic apertures smaller than the emission region will 
also produce unequal emission and absorption equivalent 
width. Hence, it is not surprising that the LRIS spec- 
tra rarely show equal absorption and emission equivalent 
widths. 

The point is that resonance emission should be ex- 
pected in galaxy spectra, and we must consider how much 
it fills in the intrinsic absorption profile, a phenomenon 
we call emission filling. In spectra of nearby galaxies, 
the near absence of emission filling in Na I AA5890, 96 
has long been attributed to the high gas densities in re- 
gions with a significant neutral sodium fraction; the large 
pathlcngth of scattered photons produces a high proba- 
bility of absorption by a dust grain (Heckman et al. 2000; 
Martin 2005, 2006; Chen etal. 2010). Quite often, the 
Mg II regions in galaxy spectra exhibit a P Cygni profile, 
which clearly shows resonance emission filling in the part 
of the intrinsic absorption trough near the systemic ve- 
locity (Martin & Bouche 2009; Weiner et al. 2009; Rubin 
etal. 2011a; Coil etal. 2011; Erb etal. 2012.). 

Absorption from the ground state necessarily leads to 
resonance emission when the ground state has a single 
level, such as in the cases of Na I, Mg II, and C IV ions, 
but can be followed by fluorescent emission when the 
ground state has fine structure, such as for the Fe + and 
Si + ions. In the composite spectra in Figure 2, we iden- 
tify fluorescent Fe II* lines at A2365.56, 2396.26, 2612.56, 
and 2626.45; the Fe II* A2632.10 line may be marginally 
detected in a few individual spectra. Because these Fe II* 
emission lines all have upper levels that are populated 
by absorption from the ground state, we can attribute 
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them to fluorescence. No absorption from excited levels 
in Fe II is detected consistent with densities well below 
the critical densities of the Fe II transitions. 

We understand the reduction in Fe II A2383 equiva- 
lent width relative to Fe II A2374 in terms of emission 
filling. The spectroscopic selection rules for total an- 
gular momentum imply that the Fe II A2383 photons 
will scatter like the common transitions (such as Mg II 
AA2796, 2803 ) with a singlet, ground state. Absorption 
of Fe II A2374 photons, however, can produce fluorescent 
Fe II* A2396.26 photons that do not fill in the resonance 
absorption trough. The relative radiative decay rates 
suggest that 88% of Fe II A2374 absorptions will decay 
as Fe II* A2396.26 photons. 

Fluorescent emission following Fe II A2344 absorption 
in Fe II* A2381.49 will also contribute to filling in the 
Fe II A2383 absorption trough, but the resonance scat- 
tering of the latter will dominate the emission filling of 
the Fe II A2383 trough. An analogous situation arises in 
the longer wavelength Fe II multiplet shown in the mid- 
dle panel of Figure 3, where absorption in the transition 
with the higher oscillator strength, Fe II A2600, usually 
results in emission of a photon at the same wavelength. 
In contrast, the excited level of the weaker transition 
Fe II A2587 typically decays via a longer wavelength, and 
therefore fluorescent, photon. Comparison of the Ein- 
stein A coefficients among all the transitions from the 
excited level indicates a decay via Fe II* A2612.65 45% 
of the time, through Fell* A2632.10 23% of the time, 
and in the resonance transition A2587 the remaining 32% 
of the time. These fluorescent emission lines, and the 
Fe II* A2626 fluorecence that follows the absorption in 
the Fe II A2600A resonance transition (13% of the time), 
are clearly detected in the spectra shown in Figures 2- 
5 and will be discussed further in a forthcoming paper 
(Kornei etal. , in prep.; Martin etal. , in prep.). 

Photons absorbed in the Mg I A2853 transition will 
scatter (rather than fluorescence) because the ground 
state is again a singlet. Could emission filling also ex- 
plain the decreasing equivalent width of Mg I A2853 in 
the spectra of lower mass galaxies seen in Section 4. It 
would be premature to draw this conclusion from the 
equivalent widths of the Mg I A2853 feature alone; any of 
a number factors - such as a decrease in the gas velocity 
dispersion, an increase in the Mg II ionization fraction, 
or a reduced covering fraction - would systematically de- 
crease the equivalent width with decreasing stellar mass. 

The presence of Mg II emission in galaxy spectra (Mar- 
tin & Bouche 2009; Weiner et al. 2009) can often be at- 
tributed to resonant scattering (Rubin etal. 2011; Coil 
etal. 2011; Erb etal. 2012). Emission that fills in part 
of the absorption trough greatly complicates the inter- 
pretation of the Doppler shift and equivalent width of 
Mg II absorption in galaxy spectra. Radiative transfer 
models, however, predict little emission at the highest 
outflow speeds (Prochaska etal. 2011). The maximum 
blueshift of the Mg II A2796 absorption is therefore the 
only property of the intrinsic absorption profile that we 
can routinely measure directly. The maximum blueshift 
of the absorbing gas has also been used previously to dis- 
cuss whether outflowing gas is bound to the galaxy (e.g., 
Heckman etal. 2000; Martin 2005; Martin & Bouche; 
Weiner et al. 2009; Coil et al. 2011; Heckman et al. 2011). 



3.3. Fitting Fe II Absorption Troughs 

Examination of our spectra showed that emission fill- 
ing complicates the direct measurement of intrinsic ab- 
sorption properties, especially in the Mg II AA2796, 2803, 
Fe II A2383, and Fe II A2600 transitions. First, the equiv- 
alent widths measured for transitions from a single ion do 
not always follow the curve of growth (Spitzer 1978). In 
the Figure 2 spectra, the highest oscillator strength Fe II 
transition at A2383 has a lower equivalent width than 
other Fe II transitions. Secondly, a partially filled ab- 
sorption trough has a centroid velocity that is bluer than 
the centroid of the intrinsic absorption profile because 
the resonance line emission is closer to the systemic ve- 
locity than is the absorption. This prevents robust mea- 
surement of the Doppler shift without detailed modeling 
(Prochaska etal. 2011). Figures 7 and 8 in the recent 
work of Erb etal. (2012) further illustrate these points 
quantitatively. 

We circumvent these problems by adopting the cen- 
troid velocity of Fe II absorption rather than Mg II as 
our primary indicator of the net Doppler shift of the low- 
ionization gas. In our discussion of emission filling, we 
explained why we expect the Fe II A2344, 2374, and 2587 
absorption troughs to reflect the outflow opacity more 
closely than does the Mg II profile. Absorption in these 
transitions usually leads to fluorescence (i.e., emission 
of a longer wavelength photon) rather than scattering, 
so the intrinsic shape of the absorption trough is much 
less affected by emission filling than the Mg II or Fe II 
A2383 lines. We therefore fit only those transitions that 
have a high probability of fluorescence, i.e., decaying via 
a longer wavelength photon rather than scattering and 
thereby filling in the absorption trough. 

We detect resonance absorption lines at high signifi- 
cance in the individual LRIS spectra. Since the spec- 
tral resolution does not fully resolve the shapes of most 
absorption troughs, we can only robustly measure the 
Doppler shifts and equivalent widths of these absorption 
troughs in general. Centroid velocities have been widely 
used to describe galactic outflows in the past, but we 
caution the reader that interstellar absorption near the 
systemic velocity blends with the absorption profile from 
outflowing (inflowing) gas in a low resolution spectrum. 
In Section 3.3.2, we will examine how interstellar ab- 
sorption at the systemic velocity biases the estimates of 
outflow properties. 

3.3.1. Single- Component Fits Describing the Absorption 
Centroid Velocity 

In Figure 4, we show examples of the Fe II absorption- 
line profiles from our LRIS data. The Doppler shift was 
jointly fit to five Fe II transitions: Fe II A 2250, 2261, 
2344, 2374, and 2587. The Fe II transitions at A2383 and 
A2600A are excluded from the fit because they show the 
strongest emission filling. Even though they significantly 
constrain the column density, the weak Fe II 2250 and 
2261 transitions provide little constraint on the Doppler 
shift. The fitted model simply describes the net Doppler 
shift of the intrinsic Fe II absorption. 

The shapes of the absorption troughs are strongly influ- 
enced by the instrumental response and are therefore of- 
ten well-described by a Gaussian line profile. We adopted 
a more physical description of the line profile in order to 
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Fig. 3. — Comparison of the Fe II and Mg II line profiles in average spectra of galaxies with different properties. Top: Towards lower 
mass galaxies, the Mg II absorption troughs become shallower, the Mg II emission becomes stronger, and the Fe II* emission in A2612, 2626 
increases. In the highest mass tertile, the Mg II doublet ratio is not inverted, and Fe II A2383 has a higher equivalent width than Fe II 
A2374, as it should in the absence of emission filling. Middle: We see similar trends with color because color is strongly correlated with 
stellar mass (redder galaxies are more massive on average). Bottom: The resonance absorption troughs become deeper in more luminous 
galaxies. See the text for a full explanation of why the stronger absorption in more luminous, redder, more massive galaxies is largely 
attributable to emission filling. 



properly model the relative line strengths in the Fe II se- 
ries. For each transition, j, our fitting function takes the 
form 



7 J -(A) = 7 cJ (A)e 
The optical depth distribution, 



-7v (A) 



T o,j 



J -(A-Ao, j ) 3 c 2 /(6 3 Ag >;f ) 



(5) 



(6) 



introduces the Doppler parameter b in order to describe 
the location of the absorbing clouds in velocity space rel- 
ative to line center. The central wavelengths of the five 
Fell transitions, Ao,j, are tied together to define a sin- 
gle Doppler shift, V\. The optical depths of the different 
transitions are tied together by the ratio of their oscil- 
lator strengths (and the rest wavelengths of the transi- 
tions), so one parameter To(A2344) is sufficient to define 
the absorption optical depth at line center for each tran- 
sition. Equivalently, the product of the column density 
and the covering factor, N(Fe + )Cj 1 , determines the op- 
tical depth in each transition. When line profiles from 
multiple transitions overlap at a given wavelength, e.g. 
rcdshifted Fe II A2374 and blueshifted Fe II A2383, we 



add their respective contributions to the total optical 
depth at that wavelength, r(A). Provided the form of 
the optical depth distribution is symmetric about the 
centroid of the line profile, Ao, the shape of this distri- 
bution does not affect the fitted central wavelength Ao, 
which defines the Doppler shift Vi. 11 
The free parameters Vi, b, and Tn(A2344) were fit- 

11 In our schematic picture, the angular extent of the stellar 
population producing the continuum emission is much larger than 
that of a gas cloud. This partial covering of the continuum source 
has often been described by a covering factor, Cf, which represents 
the fractional area of the continuum source covered by gas clouds. 
Since a fraction (1 — Cf) of the galaxy continuum is observed di- 
rectly without any absorption, a widely used, alternative parame- 
terization of the line profile is I(X)/I C (A) = l-Cf +C/e _T ( A ) (e.g., 
Rupkc et al. 2005). This parameterization of the absorption trough, 
however, does not capture the strong velocity dependence of the 
covering factor found by Martin & Bouche (2009). For this reason, 
and because velocity-dependent covering factors (Arav et al. 2005) 
introduce structure in the line profile that cannot be observed at 
the resolution of our spectra, we simply fitted the 3-parameter 
model of Eqn. 5. For any symmetric Ot(v) distribution, partial 
covering will clearly not impact the fitted Doppler shift. In our 
simple parameterization of the line profile, Eqn. 5 where Cf = 1, 
the width of the absorption trough depends on both b and tq. 
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Fig. 4. — Four comparisons of single-component (top panels) and two-component (bottom panels) fits to the Fell series with residuals 
shown at the bottom. Top Left Panels: The single-component fit to the 600 lines mm -1 spectrum of 32011099 yields an outflow velocity of 
—46 ± 14 km s — 1 . In a two-component fit, a systemic component at zero velocity describes much of the total absorption equivalent width; 
but the fit also requires a blueshifted Doppler component (Vd op = —101 ± 25 km s — 1 , Woop = 1.16A) to model the net blueshift of the 
absorption trough. Top Right Panels: The single-component fit to the 600 lines mm -1 spectrum of 22028473 yields an inflow velocity of 
42 ± 10 km s — 1 . In a two-component fit, a systemic component at zero velocity describes much of the total absorption equivalent width; 
but the fit also requires a redshifted Doppler component (Voop = 133 ± 47 km s — 1 for the inflow and Woop = 0.46A) to model the net 
redshift of the absorption trough. Bottom Left Panels: The single-component fit to the 400 lines mm -1 spectrum of 12015563 yields an 
outflow velocity of —98 ±42 km s — 1 . In a two-component fit, a systemic component at zero velocity describes much of the total absorption 
equivalent width; but the fit also requires a blueshifted Doppler component (Voop = —154 ± 84 km s — 1 , Woop — 1.55 A) to model the 
net blueshift of the absorption trough. Bottom Right Panels: The single-component fit to the 400 lines mm -1 spectrum of 12012817 yields 
Vi = 29 ±44 km s — 1 and does not require a net flow. In a two-component fit, a systemic component at zero velocity describes much of the 
total absorption equivalent width; we classify the Doppler component as insignificant because the equivalent width is extremely low and 
the change in the fit statistic is tiny. In general, the two component fits increase the estimated velocity of the flow, relative to the single 
component fit. 



ted as follows using custom software. Each galaxy 
spectrum was de-redsfiiftcd to the rest frame and 
normalized by the continuum. A model with line 
profiles centered at zero velocity was created assum- 
ing b = 200 km s _1 and to = 1. The model 
was convolved with a Gaussian profile of width a = 
120 (185) km s _1 , which describes the instrumental 
response for the 600 lines mm~ (400 lines mm - ) spec- 
tra. The residuals between the model and the data were 
minimized iteratively using the Levenberg-Marquardt al- 
gorithm (Press et al. 1992) to make successive parame- 

When we attempted to included a covering factor in the descrip- 
tion of the line profile, we found b, Cf, and to were all strongly 
covariant due to the low resolution of our data. 



ter estimations; the median reduced chi-squared value 
was close to unity. The resulting covariance matrix pro- 
vides a good estimate of the uncertainty in a fitted pa- 
rameter when the errors are normally distributed. We 
tested this assumption by bootstrap resampling typical 
spectra. The resulting distributions of Doppler shift 
measurements were found to be normally distributed, 
so we use the corresponding element of the covariance 
matrix to estimate the 68% confidence interval on V\, 

5V\ = \J <?yi yi- The Doppler shifts from the joint fit, 

V\, were also found to be consistent with the average 
of the velocities measured in the individual fits to Fe II 
A2344, A2374, and A2587. The error distributions for 



11 



To and b, in contrast, are often not normally distributed, 
so we do not use these fitted parameters in our analysis; 
we choose instead to discuss measurements of equivalent 
width and (in a later section) the bluest absorption ve- 
locity, V max . 

These joint fits to five Fe II lines describe the Dopplcr 
shift of the entire absorption trough with one velocity, so 
we will refer to them as the single- component fits. They 
avoid transitions with strong emission filling and beat 
down statistical error by combining information from 
multiple transitions. Within the LRIS sample of 208 
spectra, 36 had no good Fe II lines to fit due sky line 
residuals (5 cases) or low equivalent widths (31 spec- 
tra); and the fit failed to converge for an additional 7 
spectra. Since the 43 spectra without a fit have low 
continuum S/N ratio, the absence of significant absorp- 
tion does not provide much information; and we drop 
these from the analysis of V\ Dopplcr shifts. Among 
165 single-component fits, we find significant blueshifts, 
Vi/S(Vi) < —3, for 35 spectra and significant redshifts, 
Vi/5(Vi) > 3, for 11 spectra. 12 When the significance 
of the Fe II Doppler shift is greater than 3cr, we list the 
properties of the galaxy in Table 2. Table 3 and Ta- 
ble 4, respectively, list the fitted Doppler shifts (column 
5) for these outflows and inflows as well as the Fe II A2374 
equivalent widths (column 2) and upper limits on the 
Fe II A2261 equivalent widths (column 3) 

3.3.2. Two-Component Fits Illustrating ISM Absorption 
at the Systemic Velocity 

Having identified bulk outflows and inflows with the 
single-component fitting, we next model the contribu- 
tion of the interstellar medium (at the systemic velocity) 
to obtain an estimate of the velocity that characterizes 
the absorption equivalent width produced by the bulk 
flow. In spectra of nearby star-forming galaxies, inter- 
stellar gas is generally the primary source of absorption 
in resonance lines near the systemic velocity (Heckman 
ctal. 2001; Martin 2005; Martin 2006; Rupke etal. 2005; 
Soto etal. 2009). Stronger interstellar absorption at the 
systemic velocity is observed towards low redshift galax- 
ies with higher inclinations (Chen etal. 2010). In these 
resonance lines, the equivalent width of the interstellar 
absorption only becomes weaker than the stellar photo- 
spheric absorption in post-starburst galaxies, which can 
be recognized by the presence of absorption lines from ex- 
cited states that are populated in photospheres but not 
in diffuse gas. For our purposes here, we will use the term 
maximum interstellar component in the broadest sense to 
represent all absorption at the systemic velocity, possibly 
including some contribution from stellar photospheres. 

We modeled the Fe II series, again excluding the A2383 
and A2600 lines due to emission filling, using two veloc- 
ity components - a systemic component with no Dopplcr 
shift and a Doppler component with a fitted velocity, 
Voop- Since the absorption troughs are not well resolved, 
we found that five parameter fits allow many degenerate 
solutions. To produce meaningful Doppler velocities, we 
fixed the width of the systemic component. Because the 
systemic component represents interstellar gas, we set 
its Doppler parameter to be consistent with the mea- 

12 As discussed further in Section 6, two of the inflows are not 
significant due to systematic errors. 



sured velocity dispersion of nebular gas previously fitted 
to the [O II] doublet. These photoionized, emission re- 
gions are believed to reside in a gaseous disk; and their 
velocity disperions have been directly measured from the 
[O II] AA3726,29 doublet. The doublet is resolved in the 
831 mm -1 (red channel) spectra, and the fitted Gaus- 
sian o~ has been corrected for instrumental broadening. 
We choose b sys = ^/2a{[OII\). We then fitted the maxi- 
mum interstellar absorption to the Fe II series by varying 
the optical depth at line center, obtaining the deepest ab- 
sorption trough at v — consistent with the data. With 
the interstellar component now fully specified, we next 
fitted the Doppler shift, Doppler parameter, and central 
optical depth of the Doppler component. 

This two-component fitting approach has two criti- 
cal advantages over simply reflecting the redshifted por- 
tion of the fitted absorption trough to create a model of 
the systemic absorption, as implemented for example by 
Weiner etal. (2009). First, it does not require fully re- 
solved troughs because the model is convolved with the 
instrumental response. Second, it can be used to esti- 
mate the effect of the ISM on the inferred velocity of 
redshifted absorption (gas inflows) as well as blueshiftcd 
absorption (gas outflows). We only fit two components 
when the single-component fit indicated a significant out- 
flow or inflow. The low resolution of our spectra does not 
constrain the two-component fits well enough to use them 
to determine whether the spectrum shows a net flow in 
the first place. 

Figure 4 shows examples of the two-component fits 
relative to the single-component fits. The Vd op veloc- 
ity from the two-component fit is bluer than the single- 
component Doppler shift V\, but the difference is typi- 
cally less than the error bar 5V\. The maximum single- 
component blucshift is —330 km s _1 ; but the blucshift 
of the Dopplcr component reaches —550 km s _1 in the 
two-component fits. This trend is consistent with that 
shown in Coil etal. (2011), where the impact of the sys- 
temic component on the outflow velocity is quantified. 
While the one- and two-component fits yield consistent 
outflow velocities for many spectra, there is a subset for 
which accounting for a symmetric absorption component 
at the systemic velocity significantly increases the esti- 
mated outflow velocity. It is therefore of interest to un- 
derstand which galaxies have this stronger interstellar 
absorption component. 

The equivalent width of the systemic component rises 
strongly with increasing stellar mass. Hence modeling 
the systemic component results in fitted outflow veloci- 
ties about 200 to 300 km s -1 bluer regardless of whether 
the systemic absorption is strong or weak. Only the most 
massive galaxies in our sample, however, exhibit excep- 
tionally strong interstellar absorption with W sys > 2.3A. 
In this paper, we use the V\ values to determine whether 
a net flow is present; and, when it is, we use the values 
of Vdo P and W sys to illustrate how the interstellar ab- 
sorption at the systemic velocity may influence the con- 
clusions drawn from the single-component fits. Columns 
6 and 7 of Tables 3 and 4 compile the Dopplcr parame- 
ters derived for the systemic components and the fitted 
values of the Doppler shifts, Voop- 

3.3.3. Equivalent Widths and Ionic Column Density 
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In general, the intrinsic equivalent width of each Fe II 
line in a spectrum is determined by a combination of 
the ionic column density, N(Fe + ), and the Doppler pa- 
rameter, b, called the curve of growth (Spitzer 1978). 
The equivalent width of optically-thin transitions is sen- 
sitive to Fe II column density, but the equivalent width 
of optically-thick transitions depends primarily on the 
velocity spread of the absorbing gas. The measured 
equivalent widths do not necessarily follow this theoret- 
ical curve of growth, however, because emission filling 
reduces the intrinsic area of each absorption trough in 
a given spectrum by a different amount. Here we dis- 
cuss the measured equivalent widths of a line not affected 
much by emission filling, Fe II A2374, and the upper limit 
on the equivalent width of a weak line, Fe II A2261, de- 
tected in composite spectra. Since the Fe II 2261 transi- 
tion has an oscillator strength 12.8 times lower than that 
of Fe II 2374, the equivalent width measurement and up- 
per limit together constrain the ionic column denstiy to 
order-of-magnitude accuracy. 

How accurately the ionic column density can be esti- 
mated from a single line depends on the optical depth of 
the transition. The column density at which the absorb- 
ing gas becomes optically thick, i.e., r(A) = 1, depends 
on the absorption cross section. At line center, the opti- 
cal depth can be expressed as 

T = 1.498 x 10~ 15 XofNCfb' 1 , (7) 
where / denotes the oscillator strength of the transition, 
the wavelength of the transition is in A, the column den- 
sity has units of cm -2 , and the Doppler parameter is in 
km s _1 . To illustrate how partial covering of the con- 
tinuum source affects the inferred ionic column density, 
we include a covering factor, < Cf < 1, which de- 
scribes the fraction of the continuum light covering by 
the absorbing clouds. 

The rough equality of the equivalent widths of Fe II 
A2344, A2374, and A2587 in many of the LRIS spectra 
indicates these transitions are typically optically thick. 
The intrinsic equivalent widths of these transitions there- 
fore provide lower bounds on the ionic column density. 
Our measurement of the equivalent width for the tran- 
sition with the lowest oscillator strength, Fe II A2374, 
yields the best lower limit on N(Fe + ). In this optically- 
thin limit, the linear relationship between equivalent 
width and ionic column yields a lower limit 

N(Fe+) > 6.30 x 10 14 cm" 2 W 2 Wi{k)Cj\, (8) 

where C/.i is the covering factor in units of C/ = 1. 
The summary of the outflow and inflow properties, Ta- 
bles 3 and 4 respectively, provides these lower limits on 
N(Fe + )Cf t i in column 4 for the total absorption trough, 
in column 8 for the interstellar component at the sys- 
temic velocity, and in column 9 for the Doppler compo- 
nent. 

The Fe II A2261 absorption trough is generally much 
weaker than the Fe II A2374 line, so we can be certain 
that the Fe II A2261 transition falls on the linear part of 
the curve-of-growth. In individual spectra, we measure 
an upper limit on the equivalent width of Fe II A2261. 
We place an upper bound on the ionic column density, 
N(Fe+) < 9.06 x 10 15 cm" 2 [3a(W 226 i)}Cj^ (9) 

where the 3cr upper limit on the equivalent width has 
units of A. The summary of the outflow and inflow prop- 
erties, Tables 3 and 4 respectively, provides these upper 



limits on N(Fe + )Cf i i in column 4 for the total absorp- 
tion trough, in column 8 for the interstellar component 
at the systemic velocity, and in column 9 for the Doppler 
component. 

3.4. Highest Velocity Mg II Absorption 

Resonance emission is not expected to fill in the high 
velocity wings of the the intrinsic Mg II absorption pro- 
file (e.g., Figure 5 of Prochaska etal. 2011). If the gas 
column is low at high velocities we expect to detect Mg II 
absorption even though the Fe II absorption is too weak 
to detect for several reasons. First, consider equal col- 
umn densities of Fe + and Mg + . In the optically thin 
limit, the Fe II A2374 equivalent width would be a fac- 
tor of 20 lower than the intrinsic Mg II A2796 equivalent 
width. Next, this equivalent width ratio could be higher 
because the cosmic abundance ratio of Mg to Fe is a fac- 
tor of 1.2, and Fe is more depleted onto grains than is 
Mg (Savage & Sembach 1996). Finally, since the first 
and second ionization potentials of iron are similar to 
those of magnesium, we might expect Fe II and Mg II 
ions to be present in the same phase of the wind. 13 The 
column density ratios calculated from grids of photoion- 
ization models show, however, that the relative fraction 
of Fe and Mg in the singly- ionized state is sensitive to the 
ionization parameter. Over the range of likely ionization 
parameters near galaxies, the Fe II column density can 
be similar to or significantly less than the Mg II column 
density (Churchill etal. 2003). Hence, the relative opti- 
cal depths in Fe II A2383 and Mg II A2803 can be similar, 
but the Mg II optical depth can easily be much larger for 
a range of physical reasons. 

3.4.1. Measurement Technique 

The LRIS spectra of 186 of the 208 galaxies cover the 
Mg II doublet. In 57 spectra, no significant absorption 
trough was detected at a rest wavelength of 2796A; some 
of these spectra are characterized by strong Mg II emis- 
sion that fills the absorption trough. Among the most 
prominent Fe II outflows, i.e., the 35 with 3a outflows, 
33 spectra provide Mg II coverage. We measured the 
properties of the Mg II absorption for composite spectra 
and 104 individual spectra excluding spectra that had 
weak A2796 absorption or sky-subtraction artifacts near 
the doublet. 

Some of the individual spectra show a Mg II profile 
with a blue wing. To characterize the minimum (maxi- 
mum) velocity of the Mg II absorption, we measure the 
velocity, V max , blucward (or rcdward) of the systemic ve- 
locity at which the spectral intensity is consistent with 
the continuum level (at the la level); specifically, we 
define I(X vmax ) = 1 - 5I(X vmax ), where (57(A) is the 
error spectrum. To measure V ma x, we added a Gaus- 
sian random deviate to the original value of the intensity 
(with standard deviation matched to the error spectrum 
at each wavelength), remeasured the maximum velocity 
of the absorption trough in this modified spectrum, and 
calculated the average V max and standard deviation af- 
ter 1000 iterations. The statistical error was taken to be 

13 Their first and second ionization potentials are are IP(Mg I) 
= 7.646 eV, IP(Mg II) = 15.035 eV, IP(Fe I) = 7.870 eV, and 
IP(Fe II) = 16.18 eV. 
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the standard deviation or half a pixel width, whichever 
was larger. 

Clearly, this definition is only meaningful when the 
blue wing of the line profile is well resolved. Further- 
more, comparisons of V max must be made at similar 
spectral S/N ratio. The Mg II A2796 absorption trough 
provides our most sensitive measurement of V max and 
is not blended with another transition. When the Fe II 
lines have a significant redshift, i.e., the 9 inflow galax- 
ies discussed in Section 6, we measure the the maximum 
(inflow) velocity from the red wing of Mg II 2803; but 
this measurement can be affected by emission filling. 

The absence of blue wings on line profiles in our LRIS 
spectra does not rule out the presence of an outflow due 
to sensitivity issues caused by instrumental smoothing 
and the variation in S/N ratio among spectra. The V max 
measurements should not be compared among the entire 
sample to determine which galaxies have outflows but do 
serve two useful purposes. First, values of V max blueward 
of approximately -435 km s _1 in 400 lines mm -1 spectra 
identify a resolved line wing, as does absorption blueward 
of roughly -282 km s _1 in 600 lines mm -1 spectra. Sec- 
ond, we can compare the Mg II blue wings among galax- 
ies with different properties by measuring the V max fea- 
ture in composite spectra constructed to have the same 
S/N ratio. Previous studies have done this type of anal- 
ysis at spectral resolution higher than our LRIS spectra 
(Weiner et al. 2009; Rubin ct al. 2010b). 

3.4.2. Contamination from Stellar Absorption 

The analysis of both composite and individual pre- 
sented in this paper excludes galaxies with colors placing 
them in the green valley or red sequence or K+A spec- 
tral classification. The LRIS spectra of many of these 22 
galaxies show significant stellar absorption. The stellar 
lines become stronger with the mass of the intermedi- 
ate age stellar population, so we recognized the contam- 
ination by comparing absorption features (from excited 
states rather than resonance transitions) in spectra of 
post-starburst galaxies to our spectra. The excited ab- 
sorption features are absent or weak in our spectra of 
blue-cloud galaxies. 

This rejection affects the measurement of the maxi- 
mum absorption blueshift because the stellar line profiles 
have broad wings. We measure high V max values for sev- 
eral of the 22 rejected galaxies. Because these galaxies 
are mostly in the high mass tertile of our sample, they 
generate a positive correlation between V max and stellar 
mass when we include them. Without them, the maxi- 
mum blueshift of the Mg II AA2796, 2803 absorption in 
Figure 3 does not vary significantly with stellar mass, 
£>-band luminosity, or color. 

In principle, broad stellar line profiles could be differ- 
entiated from outflowing gas by the symmetry of the for- 
mer about the systemic velocity. In practice, however, 
the Mg II A2803 line is filled in by resonance emission. 
In Figure 5 for example, visual comparison of the max- 
imum blueshifted Mg II A2374 and the maximum red- 
shifted Mg II 2803 absorption, shows that the latter is 
often clearly affected by emission filling (e.g., in the cases 
12012777, 22013210, and 32011682) leaving a significant 
excess of absorption on the blue side. Exceptions to 
this rule include the spectra of 22004858, 22005270, and 



42014138, where the symmetry might be used to argue 
for broad, stellar wings on the line profile. A definitive 
answer requires fitting the spectra with stellar population 
models to determine the net difference in stellar ages and 
is beyond the scope of this paper. However, in general, 
we think it would be very hard to fill in the red wing 
from stellar absorption with resonance emission because 
the emission is narrow, unresolved in the LRIS spectra, 
while stellar absorption lines in single stellar population 
models have full widths at half maximum depth of 7-8A 
(or more at ages greater than 1.4 Gyr) and broad wings 
on the line profile. Hence, it seems clear that high ve- 
locity gas is required to produce the blue wings of many 
Mg II absorption profiles. 

4. COMPARISON OF Mg II AND Fe II ABSORPTION 
TROUGHS 

Composite spectra allow comparisons between subsets 
of galaxies at equal continuum S/N ratio. Before con- 
sidering individual spectra, it is useful to establish the 
overall trends between absorption properties and galaxy 
properties in composite spectra. Table 5 gives the veloc- 
ity and equivalent width measurements for the composite 
spectra discussed below. To make the comparison be- 
tween individual spectra more meaningful, we focus on 
400 lines mm -1 (600 lines mm -1 ) spectra with contin- 
uum S/N ratio greater than 10 (7.7) in this section, ex- 
cluding the 22 spectra with the strongest stellar absorp- 
tion (based on their classification as K+A, red-sequence, 
or green- valley galaxies). Table 6 summarizes the corre- 
lations between the absorption measurements for these 
individual spectra and the galaxy properties. 

4.1. Absorption Equivalent Widths 

4.1.1. Composite Spectra 

In Figure 3, we show composite spectra by tertiles in 
stellar mass, U — B color, and rest-frame i?-band lumi- 
nosity. The equivalent width of the Mg II absorption in- 
creases quite strongly with stellar mass, redder color, and 
higher £>-band luminosity. The lower Mg II equivalent 
width in the less massive, bluer, and intrinsically fainter 
galaxies appears to result largely from increased reso- 
nance emission filling in the absorption troughs. First, 
only the composite spectrum of the highest mass quar- 
tile of galaxies has a physical ratio of the A2796 to A2803 
equivalent widths; composite spectra of the lower mass 
quartiles have larger equivalent widths in the transition 
with the lower oscillator strength. Second, the composite 
spectrum of the lowest masses show a P-Cygni emission 
line in the Mg II A2796 transition. Third, towards lower 
mass and bluer color, the absorption is more suppressed 
in the red wing of Mg II A2803 than it is in the blue 
wing of Mg II A2796. Since models predict the strongest 
emission near (and redward of) the systemic velocity, 
this trend is consistent with emission filling reducing the 
equivalent width and maximum absorption velocity of 
the red wing. 

The Fe II multiplets in the left and center panels of 
Figure 3 also show the impact of increased emission fill- 
ing at lower mass, bluer color, and fainter luminosities. 
In spectra of lower mass, bluer, and lower luminosity 
galaxies, the equivalent width of Fe II A2383 decreases 
much more than that of Fe II A2374. In contrast to what 
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Fig. 5. — Comparison of the Mg II A2796 (red spectrum) and Fell A2374 (blue spectrum) line profiles in spectra with continuum 
S/N ratio greater than 10 at A2450. The black, vertical line marks the systemic velocity. Near Fe II A2374, note the location of Fe II* 
2365 at -1124 km s" 1 (blueward) and Fe II A2383 at 1048 km s" 1 (vertical blue line). Near Mg II A2796, note the Mg II 2803 line at 
769 km s _1 (vertical red line). The location of strong telluric lines that interfere with measurement of the blue wing are noted in yellow. 
In the spectra of 22028686, 22036688, and 42006915, resonance emission prominently fills in the Mg II absorption troughs near systemic 
velocity and likely affects the shape of the Mg II absorption troughs in many spectra. Despite emission filling, however, the Mg II lines 
generally have higher absorption equivalent width than the Fe II A2374 line, as can be seen from comparison of the blueshifted portions 
of their line profiles. A blue wing on the Mg II absorption trough is partially resolved in some spectra including 12008811, 12012777, and 
42014138. 



we observed, decreasing the Fe II column density would 
reduce the equivalent width ol Fe II A2374 faster than 
the Fe II A2383 equivalent width because the Fe II A2374 
transition, due to its lower oscillator strength, becomes 
optically thin when ro(2383) « 10. 

In Figure 6 and in Kornei etal. (2012), we show com- 
posite spectra for galaxies with lower and higher SFRs. 
Robust SFRs have only been measured for the AEGIS 
subset of our sample, so the S/N ratios of these com- 
posites are lower than those in Figure 3. The Mg II 
EW increases towards higher SFR. The equivalent width 
evidently increases more strongly with mass than with 



SFR, however, because the Mg II equivalent width de- 
clines with increasing specific SFR in the middle panel 
of Figure 6. 

4.1.2. Individual Spectra 

In Figure 5, we overlay the Fe II A2374 and Mg II 
A2796 absorption profiles. When the absorption troughs 
are weak, these line profiles can be indistinguishable as 
seen in the 12011364, 12017063, and 42013921 spectra. 
Among the better resolved line profiles, however, the ab- 
sorption trough of the weaker transition (Fe II A2374) is 
often deeper near the systemic velocity than the Mg II 
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Fig. 6. — Composite spectra for subsets of galaxies with measured SFR or significant bulk flows. Top: The higher SFR galaxies have 
higher equivalent widths in Mg II, Fe II A2600, and Fe II A2383. Middle: These absorption equivalent widths have a weaker dependence 
on specific SFR. The Mg II equivalent width is slightly lower in the half of the sub-sample with the highest SFRs. The inverted equivalent 
width ratio in Mg II clearly indicates this is due to the emission filling. Bottom: The subset of galaxies with the most significant Fe II 
blucshifts is compared to the subset with significant Fe II redshifts. The continuum S/N ratio is lower in the "inflow" spectrum (10 galaxies) 
than the "outflow" spectrum (35 galaxies). The wavelengths of the resonance absorption lines are marked with dotted lines for reference; 
and the Fe II* transtions are marked with dashed lines. 



A2796 trough. The prominent Mg II emission in spec- 
tra like 12012777, 22028686, 22036688, and 32011682 is 
clearly the reason why the Fe II trough is deeper. Hence, 
we argue that resonance emission is filling in the Mg II 
absorption troughs in spectra like 22004974, 42006915, 
and 42022307. To summarize, the most striking dif- 
ference between the Fe II and Mg II line profiles is the 
paucity of Mg II absorption (relative to Fe II) near the 
systemic velocity, which we attribute in large part to 
emission filling. 

Figure 7 shows the absorption equivalent widths mea- 
sured in these individual spectra over the full range in 
stellar mass, color, and luminosity. As anticipated from 
our inspection of the composite spectra, the Mg II equiv- 
alent width measurements for individual spectra are also 
positively correlated with stellar mass, color, and blue 
luminosity. Among these trends, the correlation with 
color is both the most significant (2.58 standard devia- 
tions from the null hypothesis) and strongest (Spearman 
rank correlation coefficient of r$ = 0.44). The Mg II 
absorption equivalent width is significantly more corre- 
lated with these galaxy properties than is the Fe II A2374 



equivalent width; i.e., comparison of columns 3 and 4 in 
Table 6 quantitatively confirms the impression obtained 
from the composite spectra in this regard. Even though 
the Mg II equivalent width is correlated with the Fe II 
equivalent width, the ratio M / (2796)/M / (2374) increases 
towards redder, more massive galaxies. 

In column 2 of Table 6, we show that the Doppler shift 
of the Fe II absorption has a negative correlation with 
specific SFR. Since our sign convention assigns negative 
velocities to Fe II blueshifts, the Fe II absorption is more 
blueshifted in galaxies that are undergoing a more sig- 
nificant episode of star formation. This correlation with 
specific SFR supports the expectation that the blueshifts 
mark galactic outflows. 

We might expect larger blueshifts to generate broader 
lines and therefore have larger Fe II equivalent width; but 
the measured correlation is in the opposite sense. The 
Doppler shift V\ (from single-component fitting) is pos- 
itively correlated with the Fe II A2374 equivalent width 
- a trend that may arise, for example, from a significant 
interstellar absorption component at the systemic veloc- 
ity. It is certainly more difficult to detect a blueshifted 
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Fig. 7.— Equivalent widths of the Mg II A2796 and Fe II A2374 
absorption troughs in individual, high S/N ratio spectra. Squares 
and circles represent the 400 lines mm and 600 lines mm 
spectra, respectively. The Mg II equivalent width is positively 
correlated with stellar mass, color, and i?-band luminosity; see 
Table 6. (To mitigate the impact of emission filling on the Mg II 
absorption measurement, we have integrated the area of the ab- 
sorption from -1000 to km s _1 .) 

component if the ISM absorption at the systemic velocity 
is strong. 

4.1.3. Implications of Emission Filling 

We have presented evidence for increased emission fill- 
ing in lower mass, lower luminosity, bluer galaxies. These 
trends may be attributed to dust attenuation because, on 
average at least, the more massive and redder galaxies are 



dustier than the lower mass and bluer galaxies. Scattered 
photons absorbed by dust grains will not emerge from the 
galaxy, and we observe absorption troughs which are less 
affected by resonance emission in spectra of higher mass 
galaxies. 

Dust is not required, however, to reduce emission fill- 
ing. Indeed, we find examples of good agreement be- 
tween the shape of the Mg II and Fe II absorption troughs 
(i.e., low emission filling) among individual spectra from 
galaxies with a wide range of properties. Viewing ge- 
ometry is another way to select for (or against) emis- 
sion. The line emission is weaker (and emission filling 
is reduced) in more collimated outflows, where the half- 
opening angle of the outflow cone is 45° for example, 
than a spherical outflow (Prochaska et al. 2011). Adding 
absorption from the ISM to the radiative transfer calcu- 
lation also reduces the number of resonance-line photons 
in the emergent spectrum. The increased photon trap- 
ping in the ISM gives each continuum photon absorbed 
in an Fe II resonance transition more chances to escape 
the galaxy as fluorescent emission, thereby reducing the 
resonance emission and leaving an observed absorption 
trough similar to the intrinsic absorption profile. 

4.2. Identifying Fast Outflows with Mg II Absorption 
In some spectra, such as 22013210, 22029066, and 
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32011682, the Mgll absorption is detected to 
nificantly higher blueshifts than is Fe II. At 
same time, the maximum blueshifts of the Mg II and 
Fe II absorption are consistent in many spectra even 
though W(Mg II2796) >> W(Fe II2374), e.g., 12008197, 
22004858, 22012723, and 42006781. Overall, we find no 
significant correlations between V max and mass (r s = 
—0.31; 1.8a) or color (r s = —0.21; 1.2cr) among the indi- 
vidual spectra, consistent with with the lack of any trend 
in the composite spectra constructed by these galaxy 
properties. In addition, the Fe II Doppler shift is not 
significantly correlated with V max (Mg II); and Vi(Fe II) 
has the stronger correlation with SFR and specific SFR. 
We summarize our observations as follows: (1) nearly all 
of the high S/N ratio spectra show a blue wing indicat- 
ing outflowing material; (2) spectra with a blueshifted 
Fe II centroid usually show Mg II absorption at highly 
blueshifted velocities; but (3) roughly half of the spectra 
with Mg II absorption at a large blueshift do not show a 
net blueshift of the centroid of the Fe II absorption. 

Among the high S/N ratio spectra, the equivalent 
widths of the Mg II absorption troughs tend to be larger 
when the absorption extends to higher blueshifts (r s — 
-0.45 at 2.64a in Table 6). We measured the Mg II A2796 
equivalent width by integrating the line profile from the 
systemic velocity blueward 1000 km s _1 . Defined this 
way, the equivalent width ranged from 0.25 to 1.0 times 
the total equivalent width integrated over the emission 
and absorption in both transitions. For the purpose 
of measuring the outflow equivalent width, and velocity 
range, this definition mitigates the impact of emission 
filling allowing a more direct comparison among spectra. 
For optically thick lines, the positive correlation between 
equivalent width and V ma x requires either larger spreads 
in velocity along the sightline (e.g., more clouds or gas 
acceleration) or higher gas covering fractions when ex- 
tremely blueshifted absorption is detected. 

Figure 5 compares the Mg II and Fe II line profiles in 
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spectra with high continuum S/N ratio. Most galaxies 
with 3cr blueshifts of the Fe II lines show a resolved blue 
wing on the Mg II profile. The blue wing in these spec- 
tra is typically detected in Mg II to velocities more than 
twice the fitted single-component Fe II velocity. In spec- 
tra where the V max (Mg II) and V\ (Fe II) measurements 
differ, the Mg II A2796 absorption is detected to larger 
blueshifts, as expected at low gas column densities. It 
is difficult to determine, however, whether the Fe II cen- 
troid fitting fails to identify some galaxies with outflow- 
ing gas. Some spectra, for example, have a well detected 
blue wing of Mg II even though the centroids of their 
Fe II profiles do not show a significant outlow; these in- 
clude 7 among the nine 3cr-inflow galaxies whose spectra 
will be discussed in a Section 6. In the absence of broad 
line wings from stellar absorption, measurements of this 
maximum absorption velocity, V max , may be more di- 
rectly related to the physical speeds in galactic outflows 
than the centroid of the Mg II absorption trough. 

Among the composite spectra shown in Figure 3, the 
wavelengths where the Fe II absorption troughs meet the 
continuum do not vary with stellar mass, -B-band lumi- 
nosity, color, SFR, or specific SFR. In Table 5, the most 
significant change in V max (Mg II) is with the Doppler 
shift of the Fe II absorption. A composite spectra 
built from individual spectra with significant (3<r) Fe II 
absorption blueshifts show Mg II absorption at larger 
blueshifts than does the typical spectrum with redshifted 
Fe II absorption. In the bottom panel of Figure 6, we 
show the Mg II absorption at large blueshifts in the com- 
posite spectrum constructed from the individual spec- 
tra with the most significant Fe II outflows. The com- 
posite spectrum of the 3<r outflow galaxies has an Fe II 
Doppler shift of —119 ± 6 km s _1 . This outflow com- 
posite shows Mg II absorption out to —901 ± 99 km s — 1 . 
The V max (Mg II) of the outflow composite is higher than 
that of the mass or luminosity composites, so we argue 
that (at least among star-forming galaxies) higher veloc- 
ity outflowing gas is the most likely reason for the blue 
wings on the Mg II profiles. In Kornei etal. (2012), we 
further show that V max (Mg II) becomes more blueshifted 
in composite spectra as the surface density of star for- 
mation increases. 

5. GAS OUTFLOWS 

Resonance absorption lines from Mg II and Fe II are 
the primary near-UV spectral diagnostics of galactic out- 
flows. At spectral resolutions of a 285 to 435 km s~ x , we 
can robustly measure the Doppler shifts and equivalent 
widths of the absorption troughs in nearly all spectra. 
A net blueshift of the transitions that decay primarily 
by fluorescence, Vi(Fe II), best identifies net outflows of 
low-ionization gas. The maximum blueshift of the ab- 
sorption in the strongest line, Mg II A2796, provides a 
complementary diagnostic of the outflow when spectrally 
resolved. 

To describe the demographics of outflows in Sec- 
tion 5.1, we introduce a methodology for computing the 
fraction of spectra showing blueshifted Fe II absorption 
or blue absorption wing in Mg II and examine how the 
outflow fraction varies with galaxy properties. We then 
discuss the outflow properties and how they depend on 
galaxy parameters in Section 5.2, and Section 5.3, re- 
spectively. 
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Fig. 8. — Histogram of continuum S/N ratio measured between 
2400A and 2500A (computed from the full sample of 208 galax- 
ies). Outflows/inflows are selected by the sign of the Doppler shift, 
Vi, fitted to the Fell absorption troughs. The more secure 3<r 
blueshifts/redshifts are denotes by the solid histograms, and a more 
generous lc cut denoted by the extended open histogram of the 
same color. Among the full sample of 208 spectra, the fraction of 
Fe II outflows is very sensitive to the sample's distribution of S/N 
ratios. Net flows at the 3(7 confidence level are extremely rare in 
spectra with SNR < 5; yet their fraction grows to about 58% at 
SNR 9 (36% with S/N ratio > 9). A significant fraction of high 
S/N ratio spectra present no net flow. 

5.1. Demographics of Outflow Galaxies 

Since blueshifts mark gas with a net outflow relative 
to the stellar system, we want to compute the fraction 
of galaxies with Fe II blueshifts in order to determine 
which galaxies host outflows. The challenge is that a 
blueshifted line detected in a high quality spectrum is not 
necessarily detectable in a lower quality spectrum. We 
illustrate this problem Section 5.1.1 and introduce a sta- 
tistical method for quantifying the prevalence of outflows 
in Appendix A. We discuss the results in Section 5.1.2. 

5.1.1. Fraction of Galaxies with Net Flows 

Figure 8 shows a histogram of continuum S/N ratio. 
As anticipated, at increased spectral S/N ratio, a higher 
fraction of net outflow (and inflow) galaxies are detected. 
The continuum S/N ratio is not significantly correlated 
with B-band luminosity, stellar mass, color, or SFR; 
hence spectral sensitivity rather than galaxy properties 
shape this growth in the outflow fraction. The contin- 
uum S/N ratio is anti-correlated with the apparent B 
magnitude of the galaxy (4<r deviation from the null hy- 
pothesis) , but the correlation is not very strong (correla- 
tion coefficient r$ — —0.32) for two reasons: the galaxies 
have a range of sizes, and the observations were obtained 
under a wide range of conditions. Since the sample se- 
lection was by apparent magnitude, the anti-correlation 
between S/N ratio and redshift is not very significant. 

Sensitivity to Doppler shifted absorption also de- 
pends on spectral resolution and the strength of ISM 
absorption. At the same S/N ratio (per pixel), a 
600 lines mm~ spectrum is more likely to yield a net 
Doppler shift than a 400 lines mm -1 spectrum because 
the uncertainty in the line centroid is smaller at higher 
resolution as described in Section 2.2.1. In addition, high 
equivalent width absorption at the systemic velocity will 
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reduce the Doppler shift of the fitted, single- velocity com- 
ponent. In practice, we cannot measure the out/inflow 
fraction at Doppler shifts less than \ Vi \ < 41 km s _1 due 
systematic errors described previously in Section 2.2.1. 
Focusing our analysis on these larger Doppler shifts 
largely removes any dependence on spectral resolution. 

In Appendix A, we describe a methodology for calcu- 
lating outflow fraction that takes into account the detec- 
tion biases introduced by the variations in S/N ratio and 
resolution among our spectra. Rather than arbitrarily 
picking a minimum Doppler shift to define an outflow, 
we consider threshold velocities ranging from the largest 
blucshift measured down to the minimum Doppler shift 
that we can reliably measure, roughly 41 km s _1 . We 
assign to each spectrum a probability that the Doppler 
shift of the low-ionization gas is larger than a thresh- 
old velocity. Finally, we obtain the outflow fraction for 
that threshold velocity by averaging these probabilities 
over the entire sample computing the uncetainty directly 
from the probability distribution. This definition of out- 
flow fraction is quite general and can be applied under a 
broad range of circumstances. For example, the Doppler 
shift may represent the fitted Fe II centroid velocity Vi , 
the Doppler shift of the blue wing of the Mg II absorp- 
tion trough V m ax, or the velocity derived from a more 
complex model. This methodology can also be applied 
to describe the inflow fraction if we use the net Doppler 
shift or velocity of the red wing of a line profile. 

In Figure 9, we show the outflow and inflow fractions 
computed for four samples. Roughly 5%, 20%, and 45% 
of the spectra in our full sample (black lines) show Fe II 
blueshifts faster than -200 km s _1 , -100 km s^ 1 , and 
-50 km s _1 , respectively. In contrast, less than 8% 
of the spectra have redshifts higher than 90 km s" 1 . 
Repeating the analysis with just the highest S/N ratio 
spectra (yellow lines), we obtain the same result. Simi- 
larly, with the sample restricted to either the lower reso- 
lution (400 lines mm -1 ) spectra or the higher resolution 
(600 lines mm - ) spectra, the outflow fraction remains 
consistent with the results obtained for the full sample. 
This test demonstrates that our results for the outflow 
fraction are insensitive to S/N ratio and spectral reso- 
lution. Because many more of the LRIS spectra have 
blueshifts than redshifts, the results for the outflow frac- 
tion are better constrained (than for inflow) and we de- 
scribe the outflow results in Section 5.1.2. 

5.1.2. Dependence of Outflow Fraction on Galaxy 
Properties 

We apply the computational method introduced here 
to compare the outflow fraction among galaxies with dif- 
ferent properties. We divide our sample into tertiles 
(thirds) by redshift, mass, color, luminosity, and SFR. 
To obtain insight about the importance of starbursts, we 
also normalize the SFR by the stellar mass; this specific 
SFR is the reciprocal of the timescale required to assem- 
ble the stellar mass at the current SFR. We present the 
fraction of galaxies with a blueshifted Fe II centroid ve- 
locity and then compare to the fraction of galaxies with 
significant blue absorption wings in Mg II. It is impor- 
tant to keep in mind that all the galaxies in our sam- 
ple are star-forming galaxies; the post-starburst, green 
valley, and red sequence galaxies are excluded from this 
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Fig. 9. — Fraction of spectra with a net Doppler-shift of the Fe II 
absorption lines (relative to nebular emission lines) . The solid lines 
show the fraction of blueshifted spectra with Vi < Vcrit] and the 
dashed lines illustrate the complementary quantity, the fraction of 
galaxies with Vi > Vcrit- The black, green (solid triangles), cyan 
(open pentagons), and yellow (asterisk) curves show the results of 
the calculation defined in Appendix A for the following four sub- 
samples, respectively: all spectra, all 400 lines mm - 1 spectra, all 
600 lines mm -1 spectra, and the highest quality spectra defined as 
S/N > 8.5(6.5) for 400 lines mm -1 (600 lines mm -1 ). The domi- 
nance of net blueshifts (65%) over spectra with redshifted absorp- 
tion (35%) is apparent at V cr it = 0; and the fraction of galaxies 
with Vi < V cr it (solid line) rises steadily as the threshold speed 
for an outflow decreases. The curves for the four subsamples agree 
within the the error bars showing the 68.27% confidence interval. 
This method of calculation does a good job of correcting for vari- 
ations in spectral S/N ratio and resolution allowing us to further 
examine the influence of galaxy properties on the fraction of spec- 
tra with Doppler-shifted absorption lines. 

analysis. 

The galaxy properties measured for DEEP2 galaxies 
were a primary motivation for our study. Stellar mass 
is perhaps the most fundamental property measured for 
the entire sample; and relative masses are determined 
to within a factor of two (Bundy etal. 2006). Measure- 
ments of .B-band luminosity and U — B color are de- 
scribed in Willmer et al. (2006) for the full sample. Star 
formation rates have been measured for 51 galaxies in our 
LRIS sample. In Kornei etal. (2012), we show that these 
extinction-corrected SFRs (derived from GALEX imag- 
ing of the AEGIS field) agree with the sum of the SFRs 
derived from uncorrected UV and 24 fim. photometry. 

In Figure 10, we show the fraction of spectra with V\ 
bluer than V cr it ■ The outflow fraction measured this way 
is clearly independent of stellar mass. The outflow frac- 
tion is enhanced among galaxies in the highest tertile by 
SFR (relative to the lowest tertile). The largest vari- 
ation with galaxy properties, however, is seen between 
the tertiles with the highest and lowest specific SFRs. 
Outflows faster than roughly 100 km s _1 are nearly 
absent at low specific SFR, -10 < log sSFR (yr -1 ) < 
—9.15. Among galaxies with high specific SFR, —8.8 < 
log sSFR (yr^ 1 ) < —8.2, the fraction of galaxies with 
Fell Doppler shifts (from -200 to -100 km s _1 ) is 
boosted by up to a factor of 1.8 relative to the entire 
subsample with measured SFR. 

In the LRIS sample, -B-band luminosity is more sig- 
nificantly correlated with SFR (6.1 a from the null hy- 
pothesis, correlation coefficient r$ — —0.87) than is any 
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Fig. 10. — Fraction of galaxies with Vi bluer than Vcrit- This measure of the outflow fraction shows no dependence on stellar mass (upper 
left) or color (upper middle). The outflow fraction is slightly higher among higher luminosity galaxies (middle left) and higher redshift 
galaxies (middle right). Over the velocity range —180 < Vi(kms _1 ) < —100, the outflow fraction among the galaxies in the highest SFR 
fertile (21 < SFR( Mq yr -1 ) < 98) is about a factor of four higher than the outflow fraction measured among the galaxies in the lowest SFR 
fertile (0.9 < SFR( M© yr" 1 ) < 10). Among the galaxies with the lowest specific star formation rates (—10.2 < log sSFR(yr x ) < —9.16), 
outflows appear to be strongly suppressed at all velocities relative to the average. The outflow fraction is enhanced by a factor of 1.8 over 
the average among the fertile of galaxies with the highest sSFRs (—8.82 < log sSFR(yr~ 1 ) < —8.25). The calculation of the outflow 
fraction follows the methods outlined in Appendix A. 



other galaxy property. 14 Hence, within the full LRIS 
sample, we might expect to find the outflow fraction el- 
evated among the more luminous galaxies. In Figure 10, 
the outflow fraction is indeed highest at all V cr it for the 
most luminous galaxies, but the significance of the dis- 
tinction is hardly compelling (< la). 

The high specific SFR galaxies have bluer colors on 
average (4.1 a from the null hypothesis, r$ = —0.58). 
Since we find a higher outflow fraction among galaxies 
with high specific SFR, we might expect an enhanced 
outflow fraction among bluer galaxies in the full sample. 
From Figure 10, we simply note that any enhancement 
in outflow fraction with bluer color is weaker than the 
marginal one seen with S-band luminosity. These results 
are consistent with outflow fraction depending primarily 
on specific SFR, and color being strongly dependent on 
galaxy luminosity and mass as well as specific SFR. In 
light of the increasing outflow fraction with specific SFR, 
the absence of any variation in the outflow fraction over 
the redshift range from 0.4 to 1.4 deserves further inspec- 
tion because galaxies have larger SFRs and specific SFRs 
at higher redshift. 

The three redshift tertiles for the LRIS sample with 
GALEX photometry are shown in the SFR - stellar mass 
plane in Figure 11. We show the SFR - Af* main se- 
quence (Noeske etal. 2007; Elbaz etal. 2007) fitted to 
the GALEX-detected subsample of 6102 DEEP2 galax- 
ies with stellar mass measurements. The specific SFRs of 
the galaxies with Fe II blueshifts are la larger than the 
average sSFR at the same stellar mass. Furthermore, the 
main sequence fitted to the GALEX-detected subsample 

14 A weaker correlation of somewhat lower significance is found 
between SFR and stellar mass (4.3 a from the null hypothesis, 
r s = 0.61). 



is entirely consistent with our fit to the SFRs and stel- 
lar masses estimated for the full sample of blue cloud 
galaxies in DEEP2 (Mostek etal. 2011). 15 We can there- 
fore conclude that the galaxies with Fe II blueshifts have 
slightly higher sSFR than the typical blue cloud galaxy. 

Among our subsample of galaxies with measured SFRs, 
only galaxies with sSFR > 0.8 Gyr -1 show secure, 
3(7 blueshifts. Since the slope of the main sequence is 
less than unity, the typical sSFR falls gradually with 
increasing stellar mass and this may explain why the 
Fe II blueshifts will be detected up to larger stellar 
masses at higher redshift. A sSFR threshold for Fe II 
blueshifts around 0.8 Gyr -1 would explain the paucity 
of 3a blueshifts in the low redshift tertile. Inspection of 
Figure 11, however, shows that not all galaxies with high 
sSFRs have blueshifted Fe II absorption in their spectra. 
Hence, it will be important to measure SFRs for the re- 
mainder of the LRIS sample and determine whether the 
minimum sSFR criterion remains a necessary condition 
for seeing blueshifted Fe II absorption in galaxy spectra. 

Given the distinct shapes of the Fe II and Mg II ab- 
sorption troughs in some spectra, we compare the demo- 
graphics of galaxies with blueshifted Fe II absorption and 
those with highly blueshifted Mg II absorption. Figure 12 
shows the fraction of galaxies with Mg II absorption at 

15 Mostek etal. (2011) use the AEGIS data to calibrate the re- 
sults of SED fitting to optical photometry. Comparison of this SFR 
- stellar mass relation to previous work is complicated by differences 
in both redshift range and methods used to derive stellar mass and 
SFR. For example, the Noeske et al. (2007) relation was measured 
at lower redshift 0.2 < z < 0.7, which may explain why our fit is 
about ~ 0.2 dex higher with similar slope. For our subsamples at 
0.842 < 2 < 1.349, the SFR - stellar mass relation in DEEP2 is 
remarkably shallower in slope than the Elbaz et al. (2007) relation 
fitted at 0.7 < z < 1.2. 
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Fig. 11. — Location of galaxies with Fe II blucshifts in the SFR-Af, and specific SFR - stellar mass diagram. The UV-corrected SFRs 
have been divided by 1.8 to put both the SFR and stellar mass on a Chabrier IMF. Within each redshift tertile, the bold line marks the 
main sequence of star-forming galaxies fitted to the full sample of DEEP2 galaxies with GALEX photometry; we excluded galaxies with 
U - B redder than the green valley, U - B > 1.135 - 0.032(M S + 21.18), from the fit. The slope of the best fit flattens a bit from the 
lowest redshift data (log(SFR) = 0.5419 log(Af*} - 4.7255) through the middle redshift tertile (log(SFR) = 0.4020 log(M*> - 2.9337) to 
the highest redshift spectra (log{SFR) = 0.2975 log(M») — 1.7488). The galaxies with the most robust Fe II blueshifts (solid blue squares) 
all lie on or above the main sequence. 



blucshifts larger than a threshold velocity, V C rit- Once 
again (as for Fell blueshifts), galaxies with high spe- 
cific SFR appear more likely to show highly blueshifted 
Mg II absorption in their spectra; the fraction of spectra 
with highly blueshifted Mg II absorption (as measured 
by Vmaa^Mg II)) shows no trend with redshift. In con- 
trast to the case for Fe II blueshifts, however, the outflow 
fraction derived from the Mg II blue wing shows not even 
a weak dependence on SFR alone. 

While the outflow fraction curves in Figure 12 are high- 
est for galaxies with higher masses, higher luminosity, 
and redder color, the enhancement with each of these 
properties is not statistically significant. The larger error 
bars in Figure 12 relative to Figure 10 arise from not all 
spectra showing resolved Mg II absorption troughs (and 
therefore no V max estimate). Within the measurement 
uncertainties, however, we conclude that the fraction of 
galaxies with outflows as determined by V ma x (Mg II) ap- 
pears to be less dependent on SFR than the outflow frac- 
tion determined by Vi(Fe II). Furthermore, both of these 
outflow indicators are remarkably insenstive to other fun- 



damental galaxy parameters including stellar mass, B- 
band luminosity, color, and redshift. 

5.2. Physical Properties of Outflows 

The properties of winds have generally been varied in 
cosmological simulations in order to fit observations of 
fundamental galaxy properties such as stellar content, 
SFR, and metallicity evolution (Dave, Oppenheimer, & 
Finlator 2011a; Dave, Finlator, & Oppenheimer 2011b). 
Direct measurements of outflow velocity and mass flux 
are therefore of great importance for modeling galaxy 
formation and evolution. Here we discuss the measured 
Doppler shifts and the Fe II column densities. We will 
illustrate plausible extrapolations to the wind velocity 
and mass flux but emphasize that, in order to accu- 
rately model the interplay of the various gas phases, well- 
resolved simulations of galactic outflows need to be pro- 
jected into the quantities we observe. 

The escape of metals from low-mass galaxies in the 
local universe has been established empirically (Martin 
1999; Martin etal. 2002; Tremonti etal. 2007). This en- 
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Fig. 12. — Fraction of galaxies with Mg II absorption detected at velocities larger (bluer) than a velocity V cr it- The black curve shows the 
entire sample with Mg II measurements; and the colored curves break this sample into thirds by galaxy properties. The Vmax measurements 
have not been corrected for instrumental smoothing of the spectra, so the curves are only physically meaningful at blueshifts significantly 
larger than the width of the instrumental profile which is —435 km s" 1 (vertical line in each panel) in the 400 lines mm -1 spectra. A high 
velocity blue wing is most prevalent among spectra of higher mass, higher luminosity, and redder galaxies and does not appear in spectra 
of galaxies with the lowest sSFR. 



richmcnt of the circumgalactic medium by winds prob- 
ably started very early in cosmic history based on the 
properties of metal-line systems at high-redshift (Meyer 
& York 1987; Lu 1991; Songaila k Cowie 1996; Elli- 
son etal. 2000; Schaye etal. 2003; Simcoe etal. 2006; 
Martin etal. 2010; Simcoe 2011). The main question is 
whether these outflows transport substantial mass. The 
winds clearly affect the chemical evolution of galaxies, 
but their role (if any) in creating the observed baryon 
deficit remains unclear (e.g., McGaugh etal. 2010). 

Since we use low-ionization metal lines to identify the 
blucshiftcd absorption, our survey does not probe the 
hot phase of a galactic wind. The low-ionization gas is 
entrained in the hot wind by the breakup of supershells 
(Fujita etal. 2009), the shear between the free wind and 
the galactic ISM (Heckman et al. 2000), and pre-existing, 
interstellar clouds (Cooper etal. 2008, 2009). The rela- 
tionship between this warm gas and hot wind fluid is 
critical to any extrapolation of our measurements to the 
wind properties. In the recent Hopkins etal. (2012) sim- 
ulation of a high-redshift galaxy (which has properties 
typical of our sample) however, the warm, low-ionization 
gas carries the bulk of the outflowing gas mass; hence we 
may be observing the dominant phase of the wind (by 
mass). The mass- loss rates derived from low-ionization 
gas should still be viewed as lower limits on the mass 
flux in the multiphase wind; for example, inspection of 
the other types of galaxies in Figure 4 of Hopkins et al. 
(2012) shows that the hot wind fluid can carry most of 
the mass at some outflow velocities. 

5.2.1. Outflow Solid Angle 

The measurements presented in Section 5.1.2 provide a 
statistical characterization of which galaxies have spectra 
with blucshiftcd low-ionization absorption. Among the 
various velocity measurements discussed, we have argued 



that the net Doppler shift of the Fe II lines provides the 
most robust indication of a net outflow of gas along the 
sightline. The fraction of spectra with blueshifted Fe II 
absorption decreases towards larger Doppler shifts. Pos- 
sible origins for this trend include variations in outflow 
velocity with galaxy properties, a velocity-dependent 
covering fraction, or an angular dependence in outflow 
velocity generated by the forces collimating the outflow. 

We can distinguish among these scenarios using the 
demographics of the sample. For example, some galax- 
ies with high SFR have spectra which do not show a net 
blucshift of the resonance absorption; and we will show in 
Section 5.3 that the blueshifts do not vary strongly with 
any galaxy property. Hence, we find no indication that 
the small fraction of galaxies with the largest blueshifts 
have unique physical properties and conclude that the 
outflow velocity does not vary strongly within our sample 
of blue galaxies. Furthermore, while gas covering fraction 
is likely velocity dependent as measured from resolved 
line profiles of lower redshift galaxies (Martin & Bouche 
2009), gas covering fraction does not simply explain why 
galaxies with essentially the same physical properties dif- 
fer as to whether their spectra show blucshiftcd absorp- 
tion. Over the full range in Doppler shift, the outflow 
fraction varies very little in any measured galaxy param- 
eter; the largest variations in outflow fraction are only a 
factor of two even though the stellar masses of the ob- 
served galaxies span more than 2 dex. We interpret this 
remarkable result as direct evidence for strong collima- 
tion of the ouflows at 0.4 < z < 1.4; and we predict that 
the spectra without Fe II blueshifts select galaxies viewed 
along a sightline perpendicular to the outflow direction. 

Support for this interpretation comes from the inclina- 
tions of the AEGIS galaxies. In Figure 20 of Kornei et al. 
(2012), we show that the five spectra with 3er blueshifts 
of the Fe II lines are found for galaxies viewed face-on, 
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i.e. i < 45°. In this scenario, the outflow fraction may 
be interpreted in terms of the solid angle of the outflow. 
For example, we consider a spherical outflow geometry. 
Let 9b be the half-angle of a biconical outflow, then the 
outflow subtends a solid angle 



Q = 4tt(1 - cos 6 B ), 



(10) 



and will be detected in a fraction f2/47r of the galaxies 
we observe. 

We also find a decrease in the outflow fraction as the 
Doppler shift becomes bluer, a result deserving further 
interpretation. We will simply note that if the covering 
factor of low-ionization gas decreases towards higher ve- 
locity, then the absorption troughs develop shallow blue 
wings as demonstrated by Martin & Bouche (2009) who 
had much higher resolution spectra of the Mg II doublet. 
The line profiles in the LRIS spectra lack the resolution 
required to directly measure the velocity dependence of 
partial covering. We observe a distinctly different prop- 
erty, a smaller fraction of the spectra show the centroid 
of the absorption troughs at higher outflow velocities. 
Since the troughs are just as deep as the absorption lines 
detected at lower blueshifts (in spectra of other galaxies 
with similar properties), the velocity-dependence of the 
outflow fraction in Figure 9 is not simply explained in 
terms of the gas covering fraction. Instead, we interpret 
the observation as evidence that outflows have a smaller 
opening angle at higher velocity. Smaller opening angle 
at higher velocity is qualitatively consistent with dynam- 
ical simulations of outflows where the opening angle is 
determined by the scale-height of the ISM (DeYoung & 
Hcckman 1994). 

Our study provides this first evidence for highly col- 
limated outflows at z ~ 1. In Figure 9, we show 
that the fraction of spectra with blueshifts increases 
as V cr it decreases. Specifically, we measure an out- 
flow fraction of 0.025, 0.20, 0.45, and 0.65 at V crit = 
200, 100, 50, and km s _1 , respectively. The half-angle 
of the outflow cone would be 13°, 37°, 57°, and 70°at 
velocities of 200, 100, 50, and km s _1 . Hence our re- 
sults indicate that the outflows subtend a larger (smaller) 
solid angle at slower (faster) speeds. Also, at all speeds, 
the outflow solid angle is much less than Air steradians. 
While bipolar outflows are well documented at lower red- 
shift (e.g., Heckmanetal. 1990; Chen etal. 2010; Bouche 
etal. 2012), only recently has evidence for outflow colli- 
mation emerged at redshifts from 0.09 < z < 0.9 (Bor- 
doloi etal. 2011; Kacprzak etal. 2012). 

The polar angle of 9b ~ 37° for the outflow at 
100 km s _1 derived from our results appears to be con- 
sistent with the distribution of outflowing gas in galac- 
tic halos. Bordoloi etal. (2011) mapped the Mg II ab- 
sorption at impact parameters (to background galaxies) 
within 40 kpc of star-forming galaxies at 0.5 < z < 0.9. 
Around blue, star-forming galaxies, the Mg II absorption 
equivalent widths are significantly higher at small polar 
angles, near the minor axis, than those measured closer 
to the major axis. Their results suggest the outflow cone 
makes an angle of 45° or less with the minor axis of the 
galaxy. 

5.2.2. Mass Outflow Rates 

We have nearly enough information to measure the 
mass flux in the low-ionization outflow. As described 



in Section 3.3.3, our measurements bound the Fe II col- 
umn density on both sides; this is a significant improve- 
ment over previous studies that detected outflows only 
in saturated lines, which yield only lower bounds. We 
have also shown that the largest Doppler shifts in Fe II 
are roughly 200 km s" 1 and provided some evidence 
that slower outflows may be simply viewed at higher 
inclinations (Kornei etal. 2012); this result constrains 
the outflow solid angle, ft. The ionization correction for 
Fe II is uncertain but seems likely to lie in the range 
X{Fe + ) = n(Fe + )/n(Fe) = 0.1 — 1 based on photoion- 
ization modeling (Churchill etal. 2003; Murray etal. 
2007). In the halo gas of the Milky Way, the mea- 
sured depletion of iron onto grains is much lower than 
its value in the cold disk; hence we conservatively take 
logd(Fe) = [Fe/H] = -0.69 from Table 6 of Savage & 
Sembach (1996) as our best estimate of the depiction. 
For a spherical flow launched at radius Ro, the mass flux 
is 



M = rivR N(H)m, 



(11) 



where the average mass per hydrogen atom is to = 
1.4 amu, O represents the solid angle subtended by an 
outflow with spherical geometry, v is the radial velocity 
of the ouflow, and N(H) the total hydrogen column den- 
sity. The mass flux is independent of the covering factor 
because covering factors less than unity increase the in- 
ferred column density and decrease the solid angle fl by 
the same factor. 

The galaxy spectra do not uniquely constrain the loca- 
tion of the gas along the sightline. In a radially diverging 
flow, the column density integral "down the barrel" will 
be dominated by the densest gas near the launch radius; 
hence, the linear dependence of the mass flux on the 
launch radius in Eqn. 11. Simple models relate Rq to 
the size of the starburst region (Chevalier & Clegg 1985) 
and/or several pressure scale heights in a gaseous disk 
(De Young & Heckman 1994). We therefore assume that 
Ro< 1 kpc in order to illustrate the mass fluxes. Some of 
the absorbing gas is clearly at much larger radii based on 
the detection of spatially extended Mg II emission (Rubin 
etal. 2011a; Erb etal. 2012; Martin etal. 2012b, in prep) 
and intervening Mg II absorption at impact parameters 
b < 70 kpc (Bordoloi etal. 2011) . 

Using the bounds on the Fe II column from Table 3, 
we note that the total "ISM plus outflow" column is less 
than log N(Fe + )Cf< 15.8 in half the outflow galaxies, 
and the highest upper limit is 16.27. Expressing Eqn. 11 
in terms of values consistent with the observations yields 



fn 

M = 23 Mg yr- 1 - 

Ro \ ( N(Fe+) 



v 



200 km s" 
3.16 x 10" 



1 kpc J V10 16 cm- 2 J \n{Fe)/n{H) 
0.5 \ / 0.20 



X (Fe+)J \d(Fe) 



(12) 



independent of the covering factor. For comparison, the 
median SFR in the AEGIS subsample is 12.3 M Q yr" 1 
for a Chabrier IMF (or, 21.9 M Q yr" 1 for a Salpeter 
initial mass function), indicating a mass loading fac- 
tor rj = irprji w 1.9(1.1). A lower ionization fraction 
would clearly increase our mass flux estimate. Taking 
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x{F e ) ~ 0.1, for example, raises the implied mass load- 
ing factor by a factor of 5. We will estimate lower mass 
fluxes, however, when a correction is made for the con- 
tribution of interstellar gas to the absorption equivalent 
width. 

In Section 3.3.2, we introduced a model that describes 
absorption at the systemic velocity with a maximum 
ISM component. We fit the typical properties of this 
ISM component and those of the Doppler component to 
the Fe II absorption. From Table 3, the column den- 
sity in the Doppler component is typically in the range 
14.47 < log N(Fe+)C f < 15.30, where the limits are 
the median values log Nnop(Fe + )C '/ . In the mass flux 
estimate, the lower column density of the Doppler com- 
ponent is partially offset by its larger blueshift (relative 
to the single component fit). Deprojecting the line-of- 
sight velocity based on the collimation discussed in Sec- 
tion 5.2.1, a typical outflow velocity for the Doppler 
component is Vd op = —476 km s . From Eqn. 12, 
we estimate a mass flux in the Doppler component of 
M = 2 — 11 M© yr -1 . For a Chabrier initial mass func- 
tion, we obtain r\ s=y 0.2 — 0.9, where 77 would be 1.8 times 
lower for a Salpeter initial mass function. 

Since previous mass-loss rates were derived from lower 
limits on gas column density (from saturated lines) , these 
new results show, for the first time, that the mass fluxes 
in the low-ionization gas are not orders of magnitude 
higher than the typical SFR in these z <~ 1 galaxies. 
Future work should be able to combine our constraint 
from "down-the-barrel" observations, absorption detec- 
tions towards background quasars and galaxies, and the 
extent of scattered resonance emission to improve the 
accuracy of the mass flux in the metal ions. 

5.3. Outflow Scaling Relations 

Empirical scaling relations for outflows help us to un- 
derstand which physical processes shape outflows. For 
example, for local starbursts, the nearly linear increase 
in the blueshift of Na I absorption with increasing rota- 
tion velocity (Martin 2005) motivated outflow models in 
which radiation pressure accelerates the low-ionization 
gas (Murray etal. 2005). Momentum-driven outflows, 
in general, also predict an inverse scaling between the 
mass loading parameter and galaxy velocity dispersion, 
r\ = M /SFR cx a^ 1 (Murray etal. 2005; Oppenheimer 
etal. 2010). 

The evidence presented in this paper for collimated 
outflows complicates testing the momentum-drive wind 
conjecture. The measured Doppler shifts cannot be de- 
projected into the outflow velocities without knowledge 
of the galaxy inclination. At this time, the best we can 
do is to look for trends of the median Doppler shift with 
galaxy properties. We emphasize, however, that high- 
resolution images of the sample would allow us to im- 
prove our interpretation of the spectral measurements. 

In contrast to studies of resolved galaxies at much lower 
rcdshifts, we have good statistical constraints on the halo 
masses of the galaxies in the LRIS sample. The halo 
masses are of interest because galaxy evolution models 
invoke feedback to suppress star formation in halos less 
massive than a mass floor at logM/,/ M < 11 (e.g., 
Bouche etal. 2010) and predict halo gas is mostly viri- 
alized (and therefore cold accretion and star formation 



suppressed) when log Mh/ M Q > 12, the mass ceiling. 
The observed clustering (of star- forming DEEP2 at z <~ 1 
) galaxies with Mb — 51og(/i) < —20, brighter than the 
faintest fertile of the LRIS sample, indicates a bias of 
1.28 (Coil etal. 2008). If 15-20% of those galaxies are 
satellites, the minimum halo mass is log Mh/ Mq > 11.3 
(Zheng, Coil, & Zehavi 2007), and the mean halo mass is 
log Mh/ M Q w 12.0. Roughly one-third of the galaxies 
populate halos where gas accretion should occur primar- 
ily by hot-mode accretion. Our survey is sensitive to 
changes in outflow properties below the mass floor and 
above the mass ceiling. 

5.3.1. Outflow Velocity vs. Stellar Mass 

In their study of composite spectra, Weiner etal. 
(2009) found that the Doppler shift of the most 
blueshifted Mg II absorption increased with stellar mass 
to the power 0.17. They measured the Doppler shift V\o% 
at 90% of the continuum intensity, which corresponds 
closely to our V max measurement for the median contin- 
uum S/N ratio of the LRIS spectra. The V max (Mg II) 
values measured for the LRIS sample, however, do not 
increase with M*, as illustrated in Figure 13. We found 
that V max (Mg II) does increase with M* when our anal- 
ysis includes the K+A and green valley galaxies, which 
have high stellar masses and older stars (contributing 
broad, stellar Mg II absorption). One reason our re- 
sults may differ is that our 400 lines mm" 1 spectra do 
not resolve the absorption troughs at blucshifts less than 
435 km s _1 , so we will overestimate the median velocity 
when a large fraction of spectra have blue wings just be- 
low that resolution limit. To test this idea, we identify 
the 600 lines mm -1 spectra by open squares and see that 
these V m ax values become less blueshifted with increas- 
ing M* in Figure 13. Furthermore, our measurements 
for the Fe II centroid velocity V\ , which are not affected 
much by the resolution, do not vary strongly with stellar 
mass either in Figure 13. The centroid velocities system- 
atically underestimate the outflow velocity in more mas- 
sive galaxies because the interstellar absorption at V sys is 
stronger; but when we estimated the velocity corrections 
with two-component fitting, the median blucshifts of the 
Doppler component did not reveal a significant trend 
with stellar mass. Spectral resolution is therefore not 
an obvious explanation for the discrepant result, yet it 
remains unclear whether stellar absorption significantly 
biases the Weiner etal. (2009) composite spectrum for 
high stellar mass. 

5.3.2. Do Outflows Escape? 

We follow Weiner etal. (2009) and use [Oil] velocity 
dispersions, cr<oii] (Weiner etal. 2006), to estimate the 
depths of the gravitational potentials for galaxies in the 
LRIS sample. 16 The relationship between the local es- 
cape velocity and a[oii] is sensitive to the location of the 

16 To the extent that the line-of-sight velocity dispersion mea- 
sured from the linewidth reflects the gravitational potential, the 3D 
velocity dispersion will be approximately v3 times larger. Averag- 
ing over the variation in the relation between the line-of-sight ve- 
locity dispersion, <J[oii] 1 an d the rotation velocity at different disk 
inclinations (Rix et al. 1997; Kobulnicky & Gebhardt 2000; Weiner 
etal. 2006), Weiner etal. (2009) argue that V c ~ 1.67 - 2a [OII] 
with scatter of ~ 25%. 
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Fig. 13. — Range of Doppler shifts measured vs. galaxy properties. The outflow velocities projected along our sightline show no correlation 
with stellar mass, SFR, specific SFR, B-band luminosity, or U — B color. Top Left: Maximum blueshift of Mg II A2796 absorption trough 
vs. stellar mass. Large symbols mark the median values and the ±34% range. Bottom Left: Doppler shift of Fe II absorption troughs 
vs. stellar mass. Large symbols denote the median values of the spectra with significant (> 3cr) blueshifts. Top Center: Doppler shift of 
Fe II absorption troughs vs. star formation rate. Bottom Center: Doppler shift of Fe II absorption troughs vs. B-band luminosity. Large 
symbols denote the median values of the spectra with significant (> 3cr) blueshifts. Top Right: Doppler shift of Fe II absorption troughs 
vs. specific SFR. Bottom Right: Doppler shift of Fe II absorption troughs vs. U-B color. Large symbols denote the median values of the 
spectra with significant (> 3cr) blueshifts. 



absorbing gas along the sightline to each galaxy. For a 
spherical outflow geometry (of any solid angle), much of 
the column density of an outflow viewed down-the-barrel 
is necessarily contributed by the gas near the launch ra- 
dius due to the inverse-square dilution of the gas density. 
The most conservative assumption is therefore that the 
gas lies a few pressure scaleheights above the disk at ra- 
dius R = 1 kpc. For an isothermal halo of radius 100 kpc, 
the escape velocity at R = 1 kpc is V esc = Q.la[on\- The 
median values of V esc computed this way are 281, 406, 
and 643 km s _1 for the low, middle, and high stellar 
mass tertiles; but the escape speeds are lower if the gas 
is further away. We detect absorption at V max values 
higher than these V esc estimates, see upper left panel 
of of Figure 13, and therefore support the conclusion of 
Weiner et al. (2009) that outflowing gas detected in Mg II 
may not return to the galaxy. 

The absence of observed growth in V max with stellar 
mass implies that a higher fraction of outflowing gas is re- 
cycled in more massive galaxies. In light of the significant 
uncertainties inherent in interpreting V max ~ including 
projection, sensitivity to S/N ratio, dependence on spec- 
tra resolution, and possible stellar contamination - the 
results presented here are probably not the final word on 
this scaling. In particular, we would not discount mod- 
els that predict V max oc V esc (Murray etal. 2005; Zhang 
& Thompson 2010); their results are supported by ob- 
servations of small samples of local starbursts observed 
at much higher resolution and S/N ratio (Martin 2005; 
Rupke etal. 2005). 

The Fe II Doppler shifts and the V max (Mg II) values 
suggests that outflows reach the circumgalactic medium. 
For example, we have measured Mg II absorption at 
blueshifts as high as 700 km s _1 ; and, these outflows 



would coast to 70 kpc in 100 Myr, a short enough time 
period for the host galaxy to remain blue even if the star 
formation rate declines due to the outflow. Bouche et al. 
(2006) found an anti-correlation between intervening ab- 
sorber equivalent width and host halo mass. In contrast, 
we find the largest Mg II equivalent widths in the more 
massive galaxies. Hence, the strengths of intervening 
Mg II absorption and the total Mg II equivalent width in 
galaxy spectra scale differently with galaxy mass. This 
apparent contradiction may be resolved by appealing to 
two of our results: (1) from two-component fitting (see 
Section 3.3.2), the equivalent width of ISM absorption 
(at the systemic velocity) grows with increasing stellar 
mass and (2) emission filling is stronger in lower mass 
galaxies (see Section 4.1.1). Since the resonance emis- 
sion in halo sightlines is negligible, and many halo sight- 
lines intersect little ISM, the physical affects driving the 
equivalent width - stellar mass trend in galaxy spectra 
largely do not apply to intervening absorption. 

Since our results indicate that outflows from z ~ 1 
galaxies subtend a solid angle much less than 4-7r, we ex- 
pect the intervening absorption from outflows to show a 
dependence on azimuthal angle, specifically the location 
of the sightline relative to the gas disk and outflow axis. 
Simple dynamical arguments show that the "blowout" 
will be perpendicular to the gaseous disk (De Young & 
Heckman 1994), so a distingushing property of outflows is 
that they would be detected in sightlines passing near the 
minor axis of a galaxy. New evidence for a bimodal distri- 
bution of azimuthal angles among strong Mg II systems 
(Bouche etal. 2012; Kacprzak etal. 2012) supports our 
conclusions that outflows (minor axis absorption) make 
a substantial contribution to the population of interven- 
ing absorbers. The average minor axis sightline within 
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Fig. 14. — Maximum blueshift of low-ionization absorption 
vs. SFR. Large symbols show the median values among the 
LRIS galaxies and the Weiner etal. (2009) measurement of com- 
posite spectra. The legend denotes the sample and the lines 
used to measure the maximum blueshift of the absorption trough. 
Measurements for local galaxies from Heckman etal. 2011 (aster- 
isk), Schwartz & Martin 2004 (filled triangle), Martin 2005 (filled 
square), Grimes etal. 2009 (crosses), Heckman etal. 2011 (trefoil), 
and Vazquez et al. 2004 (NGC 1705 with SFR from Marlowe etal. 
1997) . Over the relatively small range in SFR covered by our LRIS 
sample, the median V m ax(M-g II) values measured among individ- 
ual galaxies at z ~ 1 arc comparable to the maximum blucshifts 
of resonance lines in local samples of galaxies with the notable 
exception of the larger blucshifts among the galaxies with a domi- 
nant central oject (DCO). Starbursts in dwarf galaxies only show 
low-ionization absorption at significantly smaller Doppler shifts. 

60 kpc of a blue galaxy (at 0.5 < z < 0.9) shows more ab- 
sorption than a typical sightline at the equivalent impact 
parameter along the major axis (Bordoloi etal. 2011). 
This excess absorption provides the best constraint to 
date on the distance to which the Mg II clouds survive 
in the galactic outflows. 

5.3.3. Scaling with SFR, sSFR, and Mb 

In Figure 13, the Fe II centroid velocities are shown 
against SFR, sSFR, absolute B-band magnitude, and 
color. The median velocities of both the 3a outflows 
and the ler outflows show no systematic trend with any 
of these galaxy properties. We examined similar plots 
using the velocities of the Doppler component and found 
no correlations with the Doppler component. In the SFR 
and sSFR plots, the errorbars on the median are quite 
large owing to the small number of galaxies in each bin. 
This large uncertainty combined with the relatively small 
range in SFR (<~ 1 dex) could hide any underlying cor- 
relation between V\ and SFR. 

To increase the SFR baseline, we compare V max mea- 
surements to nearby dwarf galaxies and extremely lu- 
minous starbursts in Fig. 14. First, we note that our 
median V max (Mg II) measurements are consistent with 
the Weiner et al. (2009) measurements at slightly higher 
redshift. Second, when lower mass galaxies are included, 
the envelope describing the maximum outflow velocity 
increases with SFR. Excluding the galaxies with domi- 
nant central objects (DCOs), because their high veloci- 
ties may be related to AGN (Heckman etal. 2011), the 
slope is consistent with the V oc SFR a 35 relation mea- 
sured previously for local starbursts (Martin 2005). 



The true scatter in the V max — SFR scaling may be 
even larger than it appears in Fig. 14. Small Doppler 
shifts would go undetected unless the spectral resolution 
is very high, and echelle resolution has only been ob- 
tained for the nearby dwarf galaxies. Hence, the absence 
of data in the lower right corner of this plot is in part a se- 
lection effect. Moreover, projection effects only make the 
measured Doppler shifts lower than the outflow velocity, 
so the upper limit rather than the median is most in- 
teresting. The low-redshift galaxies shown in these plots 
are all starburst galaxies, and their SFR surface densities 
are typically higher than the outflow threshold (Heckman 
etal. 2003; Kornei etal. 2012). From this viewpoint, we 
argue that the absence of large Doppler shifts in dwarf 
galaxies - i.e., the paucity of points in the upper left of 
Figure 14 - is the most important feature of this scaling 
relation. Considering the multi-phase nature of winds, 
however, the apparent lack of high-velocity outflow in 
dwarf galaxies may be limited to warm, low-ionization 
gas. Based on Figure 4 of Hopkins etal. (2012), we sus- 
pect that most of the mass in a dwarf starburst outflow 
may be in the hot phase where the flow is moving at 
more than a few hundred km s _1 . If the highest ve- 
locity gas is in fact un-observable in low-ionization lines, 
then outflow velocity may not increase appreciably with 
SFR as argued, for example, based on the measured tem- 
peratures of the X-ray emitting gas (Martin 1999). 

Figure 15 compares our Vi(Fell) measurements for 
the outflow sample to other populations with starburst- 
driven outflows. Each of these samples includes spectra 
in which the resonance absorption lines show no signif- 
icant blueshift relative to the galaxy, but we plot only 
the outflow galaxies. Since the different surveys have 
different criterion for the minimum detectable Doppler 
shift, the absence of points in the lower right corner is 
not significant; but the maximum outflow speeds of each 
population likely represent the "down-the-barrel" view of 
their outflows. Taken together, these measurements show 
outflow speed increases with B-band luminosity. Since 
S-band luminosity is correlated with SFR, this result 
provides evidence for higher outflow velocities in galax- 
ies with higher SFRs. 

6. GAS INFLOWS 

In Section 1 we discused gas inflows in the context of 
a popular model whereby the gas accretion rate regu- 
lates the galactic SFR (Dekel & Birnboim 2006; Keres 
et al. 2009a; Keres et al. 2009b; Bouche et al. 2010; Dave, 
Finlator, & Oppenheimer 2012). Testing this theory re- 
quires observations of the gas inflows that reach galactic 
disks. Cosmological simulations agree that the infalling 
gas which fuels star formation is cold, in the sense that 
its temperature never exceeded 10 5 5 K, and that these 
cold flows gradually disappear at halo masses exceeding 
10 12 M due to the formation of virial shocks which in- 
crease the gas cooling time (Keres etal. 2005; Dekel & 
Birnboim 2006; Ocvirk, Pichon, & Teyssier 2008; Keres 
et al. 2009; Brooks et al. 2009; van de Voort et al. 2011a). 
The densest gas in these cold flows has a low covering 
factor but threads filaments of lower density gas with 
cross-sections much larger than a galaxy (Kimm etal. 
2010; Faucher-Giguere etal. 2011; Fumagalli etal. 2011). 
These inflowing streams carry significant angular mo- 
mentum that cause them to co-rotate with the galaxy 
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Fig. 15. — Outflow velocity of individual galaxies vs. abso- 
lute B-band magnitude. The legend identifies the samples and the 
transition(s) used to measure Doppler shift. The 3<r outflows in the 
LRIS sample have velocities intermediate to local starburst galaxies 
(Schwartz et al. 2004) and local luminous infrared galaxies (Heck- 
man etal. 2000). Local post-starburst galaxies (Tremonti etal. 
2007) and high-redshift galaxies (Pettini etal. 2001) show larger 
outflow velocities but also have significantly larger luminosities. 



and form a warped, extended disk; hence cold flows may 
be recognized kinematically by a velocity offset in the 
same direction as the rotation of the central disk (Stew- 
art etal. 2011a, 2011b). 

In contrast to the accretion rate onto halos, which is 
robustly predicted, the rate of gas accretion onto galaxies 
depends sensitively on feedback processes (van de Voort 
etal. 2011b). Winds transport supernova ejecta into the 
circumgalactic medium (Martin, Kobulnicky, & Heck- 
man 2001), so these metals mark gas that was once in a 
galaxy. The mixing between these metals and cold flows 
is challenging to resolve in numerical simulations, so the 
metallicity of accreting gas near galaxies remains quite 
uncertain. At least in some models, however, recyling of 
the metals carried to large distances by winds substan- 
tially increases the metallicity of the infalling gas (Op- 
penheimer etal. 2008; Dave etal. 2011b). In other simu- 
lations, in contrast, the cold flows are not even detectable 
in metal lines owing to their low metallicity (Kimm et al. 
2011; Fumagalli etal. 2011; Goerdt etal. 2011). 

The obvious difference between galactic outflows and 
inflows is the kinematics of the absorption-lines they im- 
print on galaxy spectra. Blueshifted (redshifted) lines in 
galaxy spectra necessarily mark outflowing (inflowing) 
gas because the absorbing material must lie on the near 
side of the galaxy. Near-UV spectra are equally sensi- 
tive to redshifted and blueshifted Fe II absorption. Be- 
cause blueshifts turn out to be more common, composite 
spectra miss the redshifted absorption, and spectra of 
individual galaxies turn out to be critical for identifying 
inflows (Sato et al. 2009; Rubin et al. 2012). Inflow detec- 
tions apparently require surveys of roughly a hundred or 
more galaxies. Both the measurements presented here 
(9 spectra out of 208) and recent work by Rubin et al. 
(2012) identify inflowing metals towards a few percent of 
z ~ 1 galaxies. 

These infalling streams may not arise from cold flows; 
the metals may mark gas recycled through galactic 



winds, tidally stripped gas, and/or material condensing 
out of hot halo gas. To provide some basis for distin- 
guishing among these processes, in Section 6.2 and Sec- 
tion 6.3 we individually examine the spectra of 3c inflow 
galaxies and estimate physical properties of the inflow- 
ing gas. In particular, we find evidence that the infalling 
gas detected in front of the galaxies marks an inflow- 
ing stream which subtends a larger solid angle than the 
galaxy on the sky. We then return to the discussion of 
the source of the inflow in Section 6.4 focusing primarily 
on the cold inflow model due to recent interest in, and 
theoretical predictions for, this particular scenario. 

6.1. Inflow Confirmation 

We identify a net inflow of low-ionization gas by a 
statistically significant redshift of the Fe II absorption 
troughs, Vi(Fell) > ia(V\). The uncertainty in this 
fitted centroid velocity does not take errors in the galaxy 
redshift determination into account. We reviewed the 
redshift determination and other sources of systematic 
error individually for the 1 1 galaxies with 3cr inflows and 
identified possible systematic errors in two of the three 
cases described here. 

All but one (22028473) of the inflow galaxies have 
strong emission lines in our LRIS spectrum; and, as 
described in Section 2.2, we measure the centroids of 
the emission lines to within ±19 km s _1 . Among 
these inflow galaxies, the average difference between this 
LRIS redshift and the DEIMOS redshift from DEEP2 
is 45 km s _1 with neither survey yielding systemat- 
ically higher (or lower) redshifts. The magnitude of 
this discrepancy is consistent with the internal error of 
41 km s _1 reported by DEEP2 from multiple observa- 
tions of galaxies and presumably arises from velocity gra- 
dients across galaxies. For the inflow galaxy 22028473, 
however, the galaxy redshift comes from the [O II] emis- 
sion in the DEEP2 DEIMOS spectrum because our LRIS 
spectrum does not cover optical, nebular lines. The fitted 
Doppler shift of the Fe II absorption, 41 ± 10 km s _1 in 
Table 4, is comparable in magnitude to the systematic 
error in the redshift determination. Although correcting 
this error is equally likely to (a) double the inflow speed 
or (b) negate the inflow detection, we drop this object 
from the inflow sample because it is not a significant 
(> 3a) detection when the error in the galaxy redshift is 
considered. 

Our adopted LRIS redshift for 22005270 exhibits one 
of the largest discrepancies from the DEEP2 redshift, 
69 km s _1 . The velocity gradient between apertures 
2 and 4 exceeds this discrepancy, so systematic differ- 
ences in slit placement may well be the cause. The H/3 
emission line in our spectrum of 22005270 appears very 
slightly redshifted, but we attribute this shift to resid- 
uals from a sky line 93 km s _1 redward of H/3 . We 
keep this object in the inflow sample. The LRIS red- 
shift accurately determines the systemic velocity for the 
purposes of comparing the velocities of the galaxy and 
the resonance absorption. Even if we adopted the ex- 
tremely conservative viewpoint that the redshift uncer- 
tainty was 60 km s -1 , the projected infall velocity of 
401 ± 43 km s^ 1 still yields a 3<r inflow detection. 

We do remove 22036194 from the inflow sample be- 
cause of the systematic error introduced in our Fe II se- 
ries fit by intervening absorption. Our 22036194 spec- 
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trum was recorded at a position angle of 25°, so the slit 
passes near the z = 0.913 galaxy 22035919 at an angu- 
lar separation of 7'.'1. Intervening Mg II absorption from 
a foreground galaxy near this redshift appears slightly 
redward of the interstellar Fe II A2344 absorption trough 
in our 22036194 spectrum. The net, blended absorp- 
tion trough is stronger and redder than the other lines 
in the Fe II series and therefore sysematically biased the 
fit. Rejection of the Fe II A2344 transition yields a fitted 
Fe II absorption velocity consistent with the interstellar 
gas in 22036194. We note, however, that 22035919 is un- 
likely the source of the intervening absorption due to the 
strength of the line and the > 50 kpc impact parameter. 
We detect prominent emission lines (near z ~ 0.831), in 
fact, from another galaxy closer to the 22036194 sightlinc 
in our 2D spectrum and suggest this object is the source 
of the intervening Mg II absorption. 

With 22028473 and 22036194 rejected, the 3cr inflow 
sample contains 9 galaxies. Figures 16 and 17 illustrate 
the 3(7 redshifts of the Fe II absorption troughs relative 
to the nebular emission. The joint fits to the Fe II series 
find a significant redshift because the absorption troughs 
of multiple transitions show a net redshift. The red- 
shifted absorption does not appear to be unique to the 
Fe II lines; the Mg II 2803 absorption trough (row 7 of 
these figures) often has a net redshift as well. Table 4 
provides the measured properties of the redshifted ab- 
sorption troughs. Table 2 summarizes the quality of the 
individual spectra and galaxy redshift, B-band luminos- 
ity, color, stellar mass, and (when available) SFR. 

6.2. Spatial Extent of Infalling Gas 

We examined the two-dimensional spectra for velocity 
gradients along each slitlit. For illustration, Figure 18 
shows the positions of the slitlets and selected regions 
of the two-dimensional spectra for the four galaxies with 
large velocity gradients. The [O II] emission lines are 
clearly tilted along the slits crossing 12011428, 12019996, 
and 42006875 as is the [O III] A5007 emission along the 
22005270 slit. We obtain additional insight about the 
inflows by investigating whether the emission or absorp- 
tion anywhere along the slit is at the same velocity as 
the redshifted Fe II absorption discovered in the galaxy 
spectra (i.e., extracted from apertures centered on the 
galaxy). We now consider the two-dimensional spectra 
of each of these galaxies in further detail. 

12011428:— Along the slit that crosses 12011428, the 
line emission is not continuous. To the north of 12011428, 
a second blob of [O II] emission is detected 1'.'14 (9.6 kpc) 
away along the slit. The image shows a hint of faint, 
extended i?-band emission in this direction. An [Oil] 
emission-line spectrum extracted at this location exhibits 
a redshift of 125 ± 15 km s _1 . The emission redshift on 
this side of the galaxy is remarkable because resonance 
absorption is seen at a similar velocity, 154 ±32 km s _1 , 
in the galaxy spectrum (extracted from an aperture cen- 
tered on the galaxy continuum). This coincidence in ve- 
locity space may indicate that the infalling stream cross- 
ing the center of the galaxy reaches the galaxy near the 
star-forming knot 9.6 kpc away. Further evidence sup- 
porting a structural connection between the redshifted 
emission and absorption is found in the grism (blue) spec- 
trum extracted on the northern side of the galaxy (aper- 



ture marked in green in Figure 18). The Fe II resonance 
absorption remains redshifted in the off-center spectrum, 
and we find no spatial gradient in the absorption veloc- 
ity across the 12011428 slit. The infalling gas apparently 
covers 12011428, and one side of the galactic disk shares 
the redshift of this infalling stream. 

22005270:— Along the slit that crosses 22005270, the 
line emission is not continuous in the two-dimensional 
spectrum shown in Figure 18. Roughly 3'.'5 southwest of 
22005270, where the slit passes near galaxy 22005066, 
we detect a spatially extended blob of [O III] emis- 
sion (l'/6 or 12.2 kpc across). Further along this slit, 
we extracted another spectrum 4'.'3 (33 kpc) southwest 
of 22005270 because the [O III] emission-line profile is 
double-peaked; the stronger, longer wavelength max- 
imum at 346 km s _1 and a weaker component at 
— 158 km s -1 . The large emission redshift near 22005066 
is remarkable because Fe II resonance absorption is seen 
at a similar velocity in the spectrum of 22005270. In 
fact, the spectrum of 22005270 shows the largest Fe II 
redshift at 401 ± 43 km s _1 . Within the uncertainties, 
this absorption velocity is marginally consistent with the 
emission component at 346 ± 19 km s _1 . Further ev- 
idence for a connection between the redshifted [O III] 
emission near 22005066 and the redshifted absorption in 
the 22005270 spectrum is apparent upon re-inspection of 
the primary spectrum towards 22005270. In both Fig- 
ure 17 and Figure 18, the bright [O III] emission in the 
22005270 spectrum falls at zero velocity by definition, 
but a weak redshifted emission component is detected at 
341 km s _1 , a redshift similar to the centroid velocity of 
the Fe II absorption. We show the extracted blue spec- 
trum where the slit passes near 22005066 in Figure 17. 
The continuum S/N ratio in this spectrum is low rela- 
tive to the blue spectrum of 22005270; but we see a net 
redshift of the Fe II and Mg II absorption. These coin- 
cidences in velocity space indicate that the infalling gas 
seen in absorption towards 22005270 is part of a physi- 
cal structure covering both 22005270 and 22005066 and 
therefore at least 33 kpc across. 

12019996:— The tilt of the [O II] emission in the 
12019996 spectrogram shows the usual signature of a 
rotating disk. The emission extends C/68 (5.6 kpc) 
to the south where the Doppler shift reaches 102 ± 
10 km s _1 . The velocity gradient between the central 
aperture and the receding side is nearly as large as the 
117 ± 25 km s _1 redshift of the Fe II absorption in the 
central aperture. The infalling gas marked by the red- 
shifted resonance absorption may therefore extend all the 
way across the 12019996 disk. 

42006875:— In the 42006875 spectrum, the tilt of the 
[O II] emission shows the usual signature of a rotating 
disk. Across the 42006875 slit, the relative redshift of the 
emission reaches 68 km s _1 l'/4 (10 kpc) to the north- 
west of the center of the galaxy. The velocity gradient 
between the center of the galaxy and the receding side is 
as large as the 51±14 km s _1 redshift of the Fe II absorp- 
tion in the galaxy spectrum. The comparable redshifts 
of the spatially extended emission and the Fe II absorp- 
tion toward the center of the galaxy suggest to us that 
the infalling gas is part of a spatially extended struc- 
ture, and this structure meets the disk roughly 10 kpc 
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Fig. 16. — Spectra with resonance absorption lines (bottom 6 rows) significantly redshifted with respect to nebular emission lines (top 
row); see also Fig. 17. The solid (red) vertical line marks the systemic velocity. Dashed vertical lines denote the wavelengths of nearby 
transitions; note the airglow line at -735 km s _1 blueward of the Mg II A2796 line in 12019996 which affects the wing of A2796 blueward of 
-359 km s — 1 . The single-component fit to the Fe II series is superposed on the Fe II A2374 profile; and the two-component fit to the Fe II 
series is plotted on the Fe II 2587 profile (with Doppler component shown by the dashed, blue line). These inflow galaxies were selected 
based on the significance of their Fe II redshift, but their spectra typically show redshifted Mg II as well. Two-component fits to the Fe II 
absorption troughs yield slightly higher inflow velocities than the one-component fits. 
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southeast of the galactic center. The Mg II absorption 
trough exhibits a red wing in the A2803 consistent with 
the infalling structure contributing to the line profile. 

How do the emission-line kinematics compare to the 
rcdshift of the Fe II absorption in the galaxy spectrum 
in the other 5 spectra? The two-dimensional spectra of 
32020468 and 12100420 do not resolve the line emission 
spatially. The absence of spatially extended emission 
is not surprising in light of their compact morphology 
in ground-based images, and 12100420 also looks com- 
pact in an HST image. Yet, as illustrated in the spectra 
shown in Figure 16, the Fe II absorption is redshifted 
123 ± 40 km s" 1 in 32020468 and 71 ± 18 in 12100420. 
In 12100420, the Mg II absorption also displays a net red- 
shift; but in 32020468, the Mg II absorption is detected 
a low significance near the systemic velocity. In the re- 
maining 3 galaxies - 12016019, 12019962, and 22004858 
- the [O II] emission along the slitlets is spatially re- 
solved and very slightly tilted. Despite the small gradient 
in the velocity of the emission, however, the spectra of 
these galaxies exhibit significantly redshifted resonance 
absorption. The Mg II absorption trough also exhibits a 
net redshift in 22004858. Hence, none of these 5 spec- 
tra show line emission at redshifts near the Fe II absorp- 
tion rcdshift. It would be valuable to obtain spectra at 
other position angles and determine whether or not these 
galaxies have any extended emission at the redshift of the 
Fe II absorption as in the case of the four objects in Fig- 
ure 18. 

To summarize, the LRIS spectra cover nebular emis- 
sion lines in all 9 inflow galaxies. We have described the 
discovery of line emission near the redshift of the Fe II 
resonance absorption in 4 of these: 12011428, 12019996, 
22005270, and 42006875. Since the velocity of the spa- 
tially offset, redshifted emission is comparable to the 
Doppler shift of the absorption measured directly to- 
wards the galaxy, the inflow may mark a spatially ex- 
tended infalling structure. To assess the probability of a 
random coincidence in velocity, we also looked for emis- 
sion features with blueshifts similar to the blueshifted ab- 
sorption. For 33 of the 35 spectra with 3a Fe II blueshifts, 
we could inspect the velocity gradient of the nebular 
emission along the slit. We found four spectra with ex- 
tended emission reaching the blueshift of the Fe II ab- 
sorption and one spectrum with emission significantly 
more blueshifted than the Fe II emission. Among the 
spectra showing outflow detections, the fraction of coin- 
cidences between emission and absorption velocities (5 
out of 33 spectra) is lower than the fraction of inflow 
spectra showing extended emission near the velocity of 
the redshifted absorption (4 matches among the 9 in- 
flow galaxies). This result supports our suggestion of a 
physical association between extended disks and satellite 
galaxies (the source of the spatially extended emission in 
the inflow spectra) and the redshifted Fe II absorption 
component. We also note that the absence of nebular 
emission (anywhere along the slit) at the velocity of the 
Fe II blueshifts further strengthens our interpretation of 
those blueshifts as evidence of outflowing gas. 

6.3. Properties of Infalling Gas 
6.3.1. Which Galaxies Have Infalling Metals? 



Nine of the 208 spectra show very significant (3cr) red- 
shifts of the Fe II absorption relative to the host galaxy. 
The low fraction of unambiguous inflows appears to be 
consistent with the 6 inflows identified recently by Rubin 
ctal. (2012) among approximately 100 objects. Here, we 
examine the properties of these 9 galaxies and their en- 
riched infall. We cannot exclude a higher inflow fraction 
at lower velocities and less significant equivalent widths; 
however, the results presented for inflows in Figure 9 in- 
dicate not more than 20% of the spectra show inflows 
with Doppler shifts larger than the systematic error of 
41 km s _1 in the determination of the systemic velocity. 

In the SFR stellar mass plane shown in Figure 11, the 
masses of 4 of these 5 inflow galaxies fall in the high mass 
tertile. Given the small number of objects, however, this 
excess is hardly significant and not supported by the full 
sample in which 4 of the 9 inflow galaxies have stellar 
masses less than logM/ M Q = 10.20, the median stellar 
mass of the LRIS sample. We have only measured the 
SFRs for 5 of the galaxies with Fe II redshifts (that lie in 
the AEGIS field). Three (two) of these 5 galaxies have 
specific SFRs lower (higher) than the average shown in 
in Figure 11, hence inflows do not depend strongly on 
specific SFR. We found 3/69, 2/68, and 4/71 galaxies 
with redshifted Fe II absorption, respectively, in the low, 
middle, and high rcdshift tcrtiles of the LRIS sample. In 
spite of the limited sample size, we therefore conclude 
that inflow galaxies may sample the full range of SFR, 
specific SFR, stellar mass, and redshift. 

The galaxies with spectra showing redshifted Fe II ab- 
sorption do not have unusual properties in any way, aside 
from having infalling gas along the sightlinc to them. 
Since the galaxies with detected inflows do not share 
some physical property, we argue that most galaxies may 
have inflows and attribute the low incidence of measured 
redshifts to the geometry of the infalling streams. Specifi- 
cally, we find 3-6% sightlines towards z ~ 1 galaxies show 
redshifted Fe II, so the covering fraction of the inflows is 
quite small. 

6.3.2. Velocity 

In Table 4, we summarize the inflow speeds fitted to 
centroids of the Fe II absorption troughs. Except for 
the extreme case of 22005270, the fitted 1-componcnt 
velocities are less than 200 km s _1 and therefore consis- 
tent with Rubin etal. (2012). The gas moves at speeds 
of material bound to galaxies. Evidence for high ve- 
locity gas can be seen in the absorption-line profiles 
shown in Figures 16, 17, and 18. The largest inflow 
velocity is detected via the Mg II absorption trough in 
22005270. In the red wings of the Fe II line profiles, ab- 
sorption at velocities exceeding 600 km s _1 is detected in 
12016019, 12019962, 12100420, 22004858, 22005270, and 
42006875. Modeling the contribution of the ISM via two- 
component fitting increases the estimated inflow velocity. 
The median inflow velocity increases from 71 km s _1 in 
the single-component fits to 191 km s _1 in the two- 
component fits. 

6.3.3. Inclination 

Rubin et al. (2012) found that 5 of their 6 inflow galax- 
ies have highly inclined (i > 55°), disk- like morphologies. 
They interpreted the high inclinations as evidence that 
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Fig. 17. — Spectra of the other galaxies with redshifted resonance absorption. Colored lines have the same meanings as defined in the 
caption to Fig. 16. In addition to the down-the-barrel spectra, i.e., a spectrum extracted from the aperture centered on a galaxy), spectra 
extracted off-center are shown for the two galaxies with spatially extended continuum emission. At locations 0'.'85 northeast of 12011428 
and 2'.'94 SW of 22005270, the Fe II Doppler shifts are 123 ± 47 and 129 ± 62 km s — 1 , respectively, consistent with the redshifts measured 
down-the-barrel (towards the centers) of those two galaxies. In contrast, the nebular emission lines, which are by definition at zero velocity 
in the down-the-barrel spectra, are prominently redshifted in the spectra extracted off center. Off the center of each galaxy, the redshift of 
the nebular emission and the Fe II absorption are very similar, consistent with a common physical location along the sightline. 
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Fig. 18. — Velocity gradients in the two-dimensional LRIS spectra of inflow galaxies. Column 1: DEEP2 (CFHT, ij-band) image with 
the l'/2 by 9'.'0 slitlet and central extraction aperture marked. The compass rose marks north and (90° counter-clockwise) east. The 
slit across 22005270 passes near a galaxy identified as 22005066 in the DEEP2 photometric catalog. Column 2: Two-dimensional LRIS 
spectra near Fe II A2600 (white tic marks). Wavelength increases from left to right, and the width of the region shown is approximately 
183 A. We denote the central apertures and the apertures used for the off-center spectra shown in Column 4 and Figure 17. Column 3: 
Two-dimensional LRIS spectra near the [O II] AA3726, 3729 doublet (or, in the case of 22005270, near [O III] A5007. Wavelength increases 
from left to right, and the width of the region shown is approximately 81 A. We mark the central apertures and the off-center apertures 
where we extracted the emission-line profiles shown in Column 5. The tilt of the emission lines reflects the rotation of dense gas in these 
galaxies; and a satellite galaxy is clearly detected in the 22005270 spectrum. Column 4 o,nd 5: Extracted spectra zoomed in on the Fe II 
A2587 and [O II] (or [O III]) lines. In all 4 spectra, the gradient in the emission velocity along the slit is as large as the redshift of the Fe II 
absorption down-the- barrel (central aperture). The infalling gas detected towards the center of each galaxy has a redshift that matches the 
Dopplcr shift of either that galaxy's outer disk or a satellite galaxy. 



rotation dominates the kinematics of the infalling gas 
but pointed out a paradox - i.e., that the absorbing gas 
must be located on only the receding side of these galax- 
ies. Kornei etal. (2012) show the morphologies of five 
of the LRIS inflow galaxies. The axis ratio was not well 
measured for 12011428. The measurements for 12100420, 
12019962, and 12019996 indicate i < 55° however, and 
only 12016019 is consistent with a large inclination of 
~ 61°. Hence, it remains unclear whether high inclina- 
tion is a necessary condition for the detection of enriched 
infall. 

6.3.4. Infalling Mass 

The (median) equivalent widths of the Fe II A2374 
and 2261 lines suggest the total Fe II column density, 
ISM plus inflow, typically lies in the range 15.08 < 



log N(Fe + )( cm~ 2 ) < 15.85. For purposes of illustra- 
tion, we conservatively choose solar metallicity gas, no 
depletion of Fe onto grains, and an ionization fraction 
of 100% singly ionized iron. The total hydrogen column 
density is then 

log N(H)( cm" 2 ) fa 
19.58 + (logiV(Fe + ) - 15.08) - log(Z/ Z ). (13) 

The implied column is log N(H) — 19.58 for solar metal- 
licity infall, but lower metallicities are plausible and raise 
the inferred hydrogen column. 

The galactic ISM clearly produces some of the Fe II 
equivalent width since our sightline directly intersects a 
galaxy. Although a true measurement of the inflowing 
column in ionized iron requires much higher spectral res- 
olution, the potential impact of the ISM on the equivalent 
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width of the inflowing component can be illustrated with 
two-component fits to the absorption trough as outlined 
in Section 3.3.2. For the Doppler component fitted to the 
Fe II absorption, we integrated the equivalent widths of 
the A2374 and A2261 profiles to obtain the lower and up- 
per bounds, respectively, on ND op (Fe + ) in Table 4. The 
median bounds provide a very rough estimate of typical 
inflow column 14.4 < log N Dop (Fe+)C f < 15.15. The 
total hydrogen column of the infalling column is therefore 
typically larger than 

log N(H)( cm" 2 ) w 

18.90 + (log N(Fe + ) — 14.40) — log(Z/ Z Q ). (14) 

Gas columns larger than those associated with galactic 
disks and damped Lyman-a systems, logA#( cm~ 2 ) w 
20.3, seem unlikely. Hence, we expect the metallicity of 
the rcdshifted gas is greater than roughly 0.04 Z Q . 

6.4. Comparison of Inflow Properties to Models 

Many properties of the infalling gas are consistent with 
the cold flow scenario. The velocity components along 
our sightlines are largely consistent with virial motion, 
the exception being the extremely large redshift towards 
22005270 which is associated with the satellite galaxy 
22005066. The two-dimensional spectra provide spatial 
information about the inflow in 3 other objects, and in 
those the redshift of the infalling gas matches the pro- 
jected velocity on one side of a rotating, gas disk. This 
association could reflect a stream crossing in front of each 
galaxy and connecting to the outer disk in the manner 
predicted by Stewart et al. (2011a, b). The low incidence 
of redshifted Fe II absorption in spectra of z ~ 1 galaxies 
would also be expectd for cold flows. 

To have any chance of detecting a cold flow in a metal 
line, the sightline must run along the filament (e.g., see 
Figure 2 in Kimm et al. 2010), and the chance of such a 
favorable orientation is small. While differences in the 
radiative transfer calculations produce some variation in 
the average covering fraction among these models, galax- 
ies with covering fractions of more than a few percent, 
as originally suggested by Dekel et al. 2009a, are now 
thought to be quite rare (Kimm etal. 2009). Numerical 
simulations suggest that only a few percent of the sight- 
lines passing within 100 kpc of a massive galaxy at z <~ 2 
intersect dense (log Nh( cm~ 2 ) > 20.3), infalling gas 
(Kimm etal. 2010; Faucher-Giguere etal. 2011; Fuma- 
galli et al. 2011). It must be acknowledged, however, that 
these high columns represent the cool gas threading the 
filaments rather than the smooth component of the cold 
streams, which dominates the cross section of neutral 
hydrogen absorption below log Nhi( cm~ 2 ) = 18 within 
Rvir (Fumagalli etal. 2011). At z < 2, less than half the 
cross section in the range 19 < log Nhi( cm~ 2 ) < 20, 
and less at higher columns, is likely due to streams (Fu- 
magalli et al. 2011). 

Theory predicts the disappearance of cold flows in ha- 
los more massive than roughly 10 12 M . The halo masses 
of many of the galaxies in the high stellar mass fertile of 
the LRIS sample arc likely above 10 12 M Q based on ci- 
ther their clustering properties (Coil et al. 2008) or halo 
abundance matching (Behroozi, Conroy, and Wechsler 
2010). In the LRIS sample, we did not find any change 
in the inflow fraction with increasing stellar mass; how- 
ever, this sample excludes red sequence galaxies. As can 



be seen by comparing the density of red sequence and 
blue cloud galaxies in Figure 1, the fraction of red galax- 
ies increases quite quickly among galaxies with stellar 
masses logM/ M Q > 10.3, the typical stellar mass in a 
10 12 M Q halo at z <~ 1 indicated by abundance match- 
ing (Behroozi etal. 2010). Whether or not the fraction 
of inflow galaxies declines with mass in an unbiased sam- 
ple cannot be determined from our data since we do not 
measure the inflow fraction towards red galaxies. Like- 
wise, the lack of any evolution with redshift might be 
limited to blue cloud galaxies which compose the LRIS 
sample. 

According to some models (Kimm etal. 2011; Fuma- 
galli et al. 2011), detection of the inflows via a metal line 
definitely implies these are not cold accretion flows. From 
our perpsective, however, the metallicity of the cold in- 
flows near galaxies remains a highly uncertain quantity 
in the models due to the unrealistic treatment of metal 
recycling via galactic winds. The ciculation of metals 
might significantly enrich the cold flows before the in- 
falling gas is incorporated into the galactic disk. 

In this recycling scenario, the distinction between cold 
flows and infalling relic outflows may become blurred. 
Numerical simulations including outflows have suggested 
that the covering factor of high column density (DLA- 
like) gas increases with halo mass, due in large part to the 
kinematics of outflowing gas (Hong etal. 2010; Faucher- 
Giguere etal. 2011). Given the lack of a strong scal- 
ing relation between outflow velocity and stellar mass, 
we expect a larger fraction of the outflowing material to 
fall back onto the more massive galaxies. Hence, the re- 
cycling of outflows should produce an increasing inflow 
column as stellar mass increases. 

The distinction between the increasing (decreasing) 
inflow column with stellar mass due to wind recycling 
(cold flows) may provide a means to discriminate the 
primary origin of infalling metal-enriched gas. This pro- 
posal further emphasizes that inflow fraction measure- 
ments should include red sequence galaxies. Although 
we expect only star-forming galaxies to drive outflows, 
the galaxies may migrate across the color - magnitude di- 
agram (due to the cessation of star formation) before the 
bulk of the outflow slows down and turns around (due to 
the gravitational attraction of the galaxy). The first in- 
flow galaxies identified appeared to support this scenario; 
their optical colors placed them on the red sequence, but 
their UV-optical colors indicated they were forming stars 
a few 100 Myr prior (Sato etal. 2009). A subsequent 
study of 13 K+A galaxies, however, discovered 2 inflow 
galaxies (Coil etal. 2011) suggesting an inflow fraction 
similar to our result for blue cloud galaxies. 17 Finally, 
mergers also bring gas into galaxies. The association of 
the galay 22005066 with the redshifted Fe II absorption 
in 22005270 may reflect this scenario. Because 22005066 
is the less luminous of the two galaxies, we would expect 
a tidal stream to pull gas out of 22005066 onto 22005270. 

In summary, cold flows, mergers, and wind recycling 
are all expected to contribute to gas infall. The low in- 
cidence of infall suggests that whichever process dom- 

17 The galaxies in our LRIS sample that show rcdshifted Fe II 
absorption do not have the spectral signatures of K+A galaxies 
nor the the colors of green-valley galaxies. They arc normal blue 
cloud galaxies and are not post-star burst galaxies. 
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inates, the covering fraction of the inf ailing streams is 
low. Much larger samples of individual spectra will be 
needed to further distinguish among these mechanisms 
on a statistical basis. 

7. SUMMARY & IMPLICATIONS 

Our understanding of galaxies will be greatly improved 
once we understood what controls the flow of matter be- 
tween galaxies and the circumgalactic medium. Empiri- 
cal constraints are essential in this regard because there 
are many physical processes that may be important and 
these processes can be nonlinear - e.g., the cooling rate 
or the interplay between mechanical and radiative feed- 
back. 

In this paper, we used unusually deep spectroscopy of 
galaxies at 0.4 < z < 1.4 to investigate which galax- 
ies have low-ionization outflows. The blue sensitivity of 
LRIS and collecting power of the Keck I telescope made 
it possible to detect NUV absorption lines in spectra of 
individual galaxies. These galaxies are typical blue cloud 
galaxies and have higher SFRs than do most galaxies to- 
day. We have fit a series of Fe II absorption lines to 
describe the bulk flow (to within ±40 km s _1 ) and con- 
strain the Fe II column density (to within an order of 
magnitude). These Doppler shifts have been compared 
to the most sensitive tracer of high-velocity gas in the 
spectra, the Mg II AA2796, 2803 doublet, which is fre- 
quently detected out to impact parameters of 70 kpc in 
galaxy halos. 

We measure Fell blueshifts of —50,-100, and — 
200 km s _1 or more in spectra of, respectively, 45%, 
20%, and 2.5% of blue-cloud galaxies at z <~ 1. This out- 
flow fraction is roughly 3 times higher in galaxies with 
21 < SFR ( Mq yr _1 ) < 98 than in galaxies with 0.9 < 
SFR ( M yr _1 ) < 10. Comparing the average spectra 
of the sample divided at the median SFR, the Fe II ab- 
sorption is significantly blueshiftcd (— 38 ± 8 km s _1 ) in 
the higher SFR, but not the lower SFR, composite spec- 
trum. Among the individual spectra, however, the corre- 
lation between Fe II Doppler shift and SFR, while in the 
direction of larger blueshifts at higher SFR (Spearman 
rank-order correlation coefficient rs = —0.225), deviates 
by only 1.6a from the null hypothesis of no correlation. 
We find the most significant outflows in spectra of indi- 
vidual galaxies with specific SFR larger than 0.8 Gyr -1 . 
Galaxies with 1.6 < sSFR ( Gyr -1 ) < 5.6 are roughly 
5 times more likely to show blueshifted Fe II absorption 
than those with 0.07 < sSFR ( Gyr -1 ) < 0.7. These re- 
sults for the outflow fraction and outflow velocity confirm 
the expectation that the blueshifts of the Fe II lines iden- 
tify outflows powered by massive stars. When combined 
with the measured cosmic SFR history, they also predict 
an increase in the fraction of galaxies with blueshifted ab- 
sorption towards higher redshift. Among just blue galax- 
ies (our LRIS sample), however, we do not detect a sta- 
tistically significant increase in the blueshifted fraction 
with redshift. 

While we find no variation of Fe II Doppler shift with 
stellar mass, blue luminosity, or U — B color, the Mg II 
equivalent width does increase towards higher stellar 
mass, i?-band luminosity, and redder color in both com- 
posite and individual spectra. Fitting the ISM absorp- 
tion explicitly shows that the increase in the Mg II equiv- 
alent width can be attributed in part to a stronger inter- 



stellar gas component in the more massive galaxies. An 
associated increase in the dust mass is consistent with 
the redder color of the more massive galaxies and reduces 
their emission in resonance lines. This reduction of emis- 
sion filling in the spectra of massive galaxies (relative to 
lower mass galaxies) amplifies the positive correlation be- 
tween the Mg II absorption equivalent width and stellar 
mass. 

We have described gas flows over a range in stellar mass 
spanning the halo mass floor and mass ceiling required 
by equilibrium SFR models. As the mass on the red 
sequence grows with cosmic time, the fraction of mas- 
sive (logM,/ M©> 10.3) galaxies having blue colors de- 
creases. It follows that although we find no change in the 
outflow properties of blue cloud galaxies between redshift 
0.4 and 1.4, the decline in the fraction of galaxies with 
blue galaxies may render outflows less prevalent among 
massive galaxies at lower redshift. 

We demonstrate that the fraction of spectra with 
blueshifted Fe II absorption is much less than 100%. 
While blueshifted absorption is found more often in spec- 
tra of galaxies with high SFRs, particularly high spe- 
cific SFRs, the blueshifted fraction is remarkably insen- 
sitive to galaxy properties overall. We conclude that the 
detection of blueshifted resonance absorption depends 
strongly on viewing angle and the outflows must sub- 
tend solid angles much less than 4-7T steradians. This 
geometry is reminiscicnt of the bipolar outflows emanat- 
ing from nearby starbursts (Heckman ctal. 1990). Be- 
cause outflows breakout of the ISM in the direction of 
least resistance, their collimation indicates the presence 
and prevalence of gas disks in blue galaxies at z <~ 1. 
The higher rotation speeds of these disks in more mas- 
sive galaxies may be the physical reason for the higher 
absorption equivalent widths observed at the systemic 
velocity in more massive galaxies. 

Based on this geometrical consideration, we can de- 
scribe the typical outflow properties. The de-projected 
outflow velocity is measured directly when our sightline 
is parallel to the outflow axis. If most of the galaxies have 
similar outflow properties, then the largest blueshifts 
we measure, roughly 200 km s _1 (line centroids) and 
500 km s _1 (Doppler components), characterize the typ- 
ical outflow speed. Using bounds on the column density 
of outflowing Fe II, the mass loss rates in the warm phase 
are comparable to the average SFR for plausible assump- 
tions about the launch radius, ionization fraction, and 
depletion. If the outflows persist at the observed rates 
for ~ 1 Gyr, then the low-ionization outfows would re- 
move roughly 10 10 M Q of warm gas, a significant bary- 
onic mass. 

The most surprising aspect of this work was the discov- 
ery of 9 galaxies with robust redshifts of low-ionization 
metals lines. Since these galaxies do not have unusual 
properties, inflows of low-ionization gas at velocities up 
to 100-200 km s _1 (based on two-component fits) appear 
to be a common property of blue galaxies at z <~ 1. We 
do not measure the hydrogen column density of these in- 
falling streams directly; but for mctallicities in the range 
of 1 to 0.1 Zq the streams would have total hydrogen 
gas cloumns of log N(H)Cf > 18.9. At these columns, 
the inflows subtend a small solid angle and escape de- 
tection in all but the largest surveys. We find some evi- 
dence that the inflows are kinematically associated with 
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satellite galaxies and extended disks, which if confirmed 
would help isolate the origin of the infalling gas. The in- 
flows may originate, for example, from tidal forces gener- 
ated by galaxy interactions or multiphase cooling of hot 
halo gas. We emphasize that aside from their metallicity, 
the properties of these inflows are consistent with recent 
predictions for cold flows of accreting gas. Since future 
modeling of the mixing between primordial inflow and 
recycled wind material may well increase the expected 
metallicity of cold flows within a few tens of kpc of a 
galaxy, the detection of the infall in metal lines does not 
obviously rule out a cold flow origin 
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APPENDIX 

OUTFLOW FRACTION CALCULATION 

As an example of how to calculate the outflow fraction, or more precisely the fraction of spectra with blueshiftcd 
resonance absorption, from data, we will use the results of our single-component fitting to the Fe II series of absorption 
lines. For each spectrum, the fitted Doppler shift, V\, and its uncertainty, SVi, provide the best estimate of the 
distribution of Doppler shifts (i.e., its mean and standard deviation) which would be obtained after many observations. 
Equivalcntly, in the special case of a symmetric probability distribution, we can describe the probability of measuring 
a Doppler shift V\ when the true velocity is equal to or past V C Ht- 

Let Pi(V < Vera) represent the probability that the absorption lines in a particular spectrum have a Doppler 
shift bluer than V cr i t . The probability, Pj(V < V cr it), is the fractional area of this distribution with Doppler shifts 
V < V cr it- Since this distribution is well described by a Gaussian function of width a = SV\ (Section 3.1), we simply 
compute Pi(V < V cr it) from an expression closely related to the normal error integral. 

The expression for the entire integral depends on the relation of V cr it and V. We adopt a sign convention that 
outflow velocities are negative. When a spectrum has a low Doppler shift relative to the outflow threshold, V cr it < V\, 
the chance that the true value of the Doppler shift is faster than V cr i t is small and given by 

Pi(V< V crit \V 1 ,dV 1 ) = 0.5-Q. (Al) 

When the fitted Doppler shift is bluer than the outflow threshold, V cr it > Vi, the chance that the true Doppler shift 
is faster than V cr u is much more significant; and the total integral is 

P{V< V crit \V 1 ,dV 1 )=0.5 + Q. (A2) 

The function 

Q{t) = -j=j\-*'*dz=\erf{t). (A3) 
describes the integral over a portion of a Gaussian distribution, where we have defined the dimensionless variable t as 

t^lVcrit-VMSVL (A4) 

For each discrete value of the threshold velocity, we compute the outflow fraction by adding the outflow probabilities 
Pi(V < V cr it) of all the spectra. We normalize this sum by the number of spectra, N spec to obtain the outflow fraction 

Pout(V crit ) = N- p l c ^rPi(V < Vcrit). (A5) 
The probability of finding a Doppler shift redder than V cr it is 1 — Pout- For positive values of V cr it, the fraction of 
spectra with V > V cr it obtained from 1 — P out describes the frequency of redshifted velocities in the sample. 

Error bars are calculated directly from the probability distribution. For a sample of N spec objects, we compute the 
probability of finding 0, 1, 2, N spe c — 1, N spec objects with V < V cr it- In practice, our code uses a subroutine that 
starts with the first spectrum in the list, calculates the probability of finding or 1 objects with V\ < V cr u, and 
then proceeds to calculate probabilities for larger samples in an iterative fashion advancing through lists of 2,3,4, ... 
objects until the probabilities are obtained for the full list of N spec objects. The sum of the first i + 1 values in this 
sequence yields the cumulative probability that the number of outflows with V\ < V cr u is less than or equal to Ni, 
where i = to N spec . The boundaries of the 68.27% confidence interval arc found by interpolating between the Ni 
values to estimate the number of spectra where the cumulative probability that V\ < V cr it is 0.1587 (lower bound) 
and 0.8413 (upper bound). 
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TABLE 1 
Kbck/LRIS Observations 



Field 


R.A. 


Decl. 


Mask PA 


Grism / Grating 


Dichroic 


Conditions 


Exposure Time 


Number of 


Observing Run 




(J2000.0) 


(J2000.0) 


(°) 


(1 mm~ 1 ) a 




Cloud Cover & Seeing(" ) 


LPJS-B/LPJS-R (s) b 


Exposures (B/R) 


Dates 


msc42_l 


02 28 56.8 


+00 33 39 


120.0 


400 / 831 


d680 


Clear to Mostly Cloudy & 0.9-1.6 


11340 / 11100 


6/6 


2007 Oct 6-7 


msc42_5 


02 30 35.6 


+00 28 17 


143.0 


400 / 831 


d680 


Mostly Cloudy to Light Cirrus & 0.7-0.9 


26883 / 25400 


15/14 


2008 Sep 28-29 


mscl2_d 


14 17 09.8 


+52 31 18 


165.0 


400 / 831 


d680 


Clear to Mostly Cloudy & 0.6-2.0 


18500 / 18000 


10/10 


2008 June 5-6 


mscl2_8 


14 18 19.4 


+52 34 33 


5.0 


400 / 831 


d680 


Clear to Partly Cloudy & 0.6-2.0 


28800 / 26940 


15/31 


2009 June 18-21 c 


mscl2_ee 


14 18 46.8 


+52 38 59 


170.0 


600 / 600 


d560 


Clear to Partly Cloudy & 0.6-2.0 


11520 / 10800 


6/12 


2009 June 18-21 c 


msc22_bb 


16 51 09.2 


+34 58 15 


25.0 


600 / 600 


d560 


Clear to Mostly Cloudy & 0.6-2.0 


18400 / 17950 


9/9 


2008 June 5-6 


msc22_6 


16 51 28.8 


+34 48 09 


60.0 


400 / 831 


d680 


Clear to Partly Cloudy & 0.6-2.0 


30760 / 23100 


16/26 


2009 June 18-21 c 


msc32_aa 


23 29 19.5 


+00 07 38 


82.0 


600 / 600 


d560 


Mostly Cloudy to Light Cirrus & 0.7-0.9 


11070 / 10800 


6/6 


2008 Sep 28-29 


msc32_l 


23 30 18.9 


+00 12 53 


80.0 


400 / 831 


d560 


Clear to Light Cirrus to Mostly Cloudy 


25810 / 24182 


14/14 


2007 Oct 6-7 & 2008 Sep 28-29 



a After binning by 2 pixels in the dispersion direction, 400 lines mm 1 and 600 lines mm -1 blue spectra have 2.18 A pix 1 and 1.26 A pix 1 . respectively. The unbinned 831 mm -1 and 600 mm 1 red spectra have dispersions of 0.93 A pix 1 
and 1.28 A pix - 1 , respectively, for all runs except 2009 June. 
^Thc typical exposure time per frame was 1800 s. 

C Thc red spectra from 2009 June were obtained with a new detector. After binning by 2 pixels in the dispersion direction, the 831 mm -1 and 600 mm -1 red spectra have dispersions of 1.16 A pix -1 and 1.60 A pix -1 , respectively. Shorter 
exposure times were used for individual frames due to high cosmic ray rates. 
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TABLE 2 

Properties of Galaxies with Net Fe II Doppler Shifts 



W DJtiCl 


Grism 


SNR Rcdshift 




T\4 n ^ ,r\cy h 




log 


or 1 V 




(1 IT) IT) \ 


Source 




1 A R tn ad I 






I ivl J 1 


(1) 
V- 1 -/ 


(2) 


(3) 


(4) 




(6) 


(7) 


(8) 


(9) 


12008197 


400 


12.2 


L 


0.9802 


-21.19 


0.804 


10.63+ 


98 


12008550 


400 


9.86 


L 


1.3024 


-21.22 


0.620 


10.08+ 


25 


12011836 


400 


9.38 


L 


0.9270 


-19.90 


0.579 


10.26+ 


14 


12012777 


400 


11.68 


L 


1.2742 


-21.06 


0.516 


10.11+ 


23 


12013242 


400 


8.83 


L 


1.2867 


-21.31 


0.450 


10.19 


26 


12011428* 


400 


4.6 


L 


1.2840 


-19.87 


0.180 


9.65 


9 


12015177 


600 


5.3 


L 


0.9860 


-20.23 


0.487 


9.69 


12 


12016019* 


400 


11.57 


L 


1.0846 


-20.71 


0.622 


10.03+ 


19 




uuu 


3.73 


£) 


1 97£zl 


91 fiR 


U. 1 uo 


10 40+ 






a nn 
4UU 


1 n 1 
1U. 1 


T 


U.0444 


i n 99 


n a c;q 


fl 77t 

y. it 1 


z 


izuiyyyo 


4UU 


7 Qft 
i .OK) 


T 

Li 


1 9Q1 9 

1 .ZolZ 


91 71 
-Zl . / 1 


n ria 


i n fi7f 


A ^ 


tot nn /i on* 


400 


fi fin 


T 

Lj 


1.1995 


-20.63 


0.601 


10.01 1 


10 




a nn 
4UU 


15.00 


L 


l.ZOoO 


91 
-Zl .04 


n A9/1 


i n x.a 
1U.04 




2200^21 (S 


400 


9.27 


L 


0.9130 


-19.49 


0.556 


9.59 




22005270* 


400 


10.57 


L 


0.8308 


-20.38 


0.710 


10.40 




22006172 


400 


10.65 


L 


0.8328 


-21.10 


0.707 


10.65 




22006207 


400 


9.42 


L 


1.2709 


-20.59 


0.419 


9.88 




22029066 


600 


10.45 


L 


0.7852 


-21.59 


0.891 


11.29+ 




22029224 


600 


8.66 


L 


0.8603 


-20.90 


0.611 


10.27 




22036854 


600 


7.29 


L 


0.9373 


-20.85 


0.564 


10.32 




22036912 


600 


6.66 


L 


0.8099 


-21.02 


0.864 


10.86 




32010773 


600 


9.23 


L 


0.8039 


-19.91 


0.234 


9.60+ 




32011098 


600 


5.01 


L 


0.9561 


-20.22 


0.546 


9.77+ 




32011099 


600 


6.51 


L 


0.8832 


-20.06 


0.492 


9.80 




32011192 


600 


9.51 


L 


0.8478 


-20.60 


0.587 


10.07+ 




32011682 


600 


8.79 


L 


0.8359 


-20.36 


0.511 


10.23+ 




32016857 


600 


8.32 


L 


0.9391 


-20.81 


0.440 


9.81+ 




32017112 


600 


5.85 


L 


1.0084 


-20.89 


0.665 


10.23+ 




32020468* 


400 


3.16 


L 


1.2355 


-21.00 


0.734 


10.28+ 




32022156 


600 


4.60 


L 


1.0434 


-20.70 


0.580 


10.10+ 




42006781 


400 


10.98 


D 


1.2859 


-20.74 


0.546 


9.68+ 




42006915 


400 


25.09 


L 


0.8945 


-20.55 


0.386 


10.50+ 




42014101 


400 


8.39 


L 


0.7494 


-20.47 


0.899 


10.94+ 




42014154 


400 


9.75 


L 


0.8427 


-19.40 


0.465 


9.77+ 




42014585 


400 


8.17 


D 


1.2705 


-20.91 


0.588 


10.18 




42014618 


400 


17.64 


L 


1.0130 


-20.37 


0.592 


9.81 




42014718 


400 


9.21 


L 


1.1904 


-21.27 


0.692 


10.85+ 




42014732 


400 


16.01 


L 


0.7502 


-19.40 


0.453 


9.50+ 




42021266 


400 


11.66 


L 


0.9765 


-19.63 


0.616 


9.50 




42022173 


400 


8.15 


D 


1.3112 


-20.25 


0.324 


9.68 




42022307 


400 


20.72 


D 


1.2595 


-21.43 


0.466 


10.22 




42026243 


400 


3.30 


L 


1.3456 


-21.66 


0.667 


10.88+ 




42033991 


400 


6.48 


L 


0.8714 


-20.13 


0.392 


9.94+ 




42006875* 


400 


19.18 


L 


0.8696 


-21.10 


0.602 


10.50+ 





Note. — (1) - Object identification from DEEP2. The Fc II absorption in these spectra is Doppler 
shifted (relative to the galaxy) at a significance level greater than 3cr. Inflows (redshifts) are denoted 
by *; the objects 22028473 and 22036194 have been dropped per the discussion in § 5.1. (2) - Blue 
grism where 400 lines mm and 600 lines mm" - 1 imply FWHM of 435 km s" 1 and 282 km s" 1 , 
respectively. (3) - Continuum SNR between 2400 and 2500 A. (4) - Source of redshift. The flag 
L/D indicates (L)RIS (the spectra presented here) or (D)EIMOS, DEEP2 redshift, respectively (5) 
- Redshift. (6) - DEEP2 absolute magnitude (Willmer etal. 2006). (7) - DEEP2 color (Willmer etal. 
2006). (8) - Stellar mass from SED fitting (Bundy et al. 2006). Galaxies with an SED fit that includes 
K-band photometry (most certain) are denoted by +; the other fits are based on BRI photometry. 
(9) - Star formation rate derived from dust-corrected UV continuum luminosity and converted to a 
Chabrier initial mass function. 



TABLE 3 
Outflow Properties 



Object 
(1) 


W(2374) 
(A) 
(2) 


W(2261) 
(A) 
(3) 


log[N(Fe+)C f ] 
( cm- 2 ) 
(4) 


Vi 

( km s- 1 ) 
(5) 


biSM 
( km s- 1 ) 
(6) 


Vdop 
( km s- 1 ) 
(7) 


log[Ar JSM (Fe+)C / ] 
( cm- 2 ) 
(8) 


log[N DO p(Fe+)C f ] 
( cm- 2 ) 
(9) 


Vmax (A2374) 
( km s -1 ) 
(10) 


Vmax (A2796) 
( km s- 1 ) 
(11) 


1 onnoi c\t 
IzUUoiy / 


l.zl ± U.zo 


0.21 


14.89 


- 15.75 


—62 ± 19 


126 


—93 ± 43 


14.80 - 


15.34 


14.58 - 


15.93 




a a.a _L i on 
— 4bD ± lzU 


1ZUU855U 


i ncr I n 1 n 

1.05 ± U.19 


0.17 


14.83 


- 15.66 


— 145 ± 23 


77 


■ ) i 1 ( 1 1 IOC 

— 322 ± 135 


14.69 - 


15.27 


14.32 - 


15.79 


—547 ± 9o 




IZUlloou 


O 1 A _l_ n QO 

Z. 14 ± U.oZ 


0.27 


15.14 


- 15.86 


— 128 ± 21 


114 


AdC _1_ 1 A O 

— 400 ± 14Z 


15.00 - 


15.65 


14.42 - 


15.01 


aciT _i_ 117 
— by / ± 117 


— do / ± y4 


1ZU1Z ( / ( 


Z. lo ± u. iy 


n 1 q 
U. lo 


15.13 


- 15.68 


— Zl 1 ± lo 


yo 


A A Q _l_ 7(1 

— 44o zt ID 


15.07 - 


15.89 


14.79 - 


15.59 


'JCxA J_ 1 OO 

— ( U4 ± 1ZZ 


nnn i 1 a p 

— yuy ± 140 


!ZUloZ4Z 


1 1/1 _L n OQ 

1.14 ± U.zo 


0.24 


14.86 


- 15.81 


— 73 ± 22 


121 


—97 ± 45 


14.89 - 


15.48 


14.78 - 


15.81 






1 On 1 K1 77 

IzUlol ( i 


O Q A _!_ n Q O 

Z.o4 ± U.oZ 


0.33 


15.17 


- 15.94 


— 101 ± 18 


91 


1 c\a _l_ a o 
— 106 ± 49 


14.95 - 


15.64 


14.72 - 


15.77 


—418 ± 51 




12019542 


0.72 ± 0.30 


0.39 


14.66 


- 16.02 


— 183 ± 60 4 


69 


—311 ± 273 


14.51 — 


15.04 


14.51 — 


15.04 






22005216 


0.98 ± 0.26 


0.25 


14.80 


- 15.82 


— 106 ± 28 


59 


—450 ± 146 


14.90 — 


15.76 


14.20 — 


14.66 




—504 d= 56 


22006172 


1.57 ± 0.28 


0.27 


15.00 


- 15.86 


— 112 ± 21 


118 


—95 ± 47 


14.84 — 


15.41 


14.80 - 


15.80 


—588 ± 178 


— 514 ± 128 


22006207 


0.76 ± 0.15 


0.16 


14.69 


- 15.63 


— 199 ± 39 


76 


—34 ± 35 


14.57 — 


15.11 


14.68 — 


16.06 


—487 ± 84 


—693 ± 94 


22029066 


1.85 ± 0.20 


0.27 


15.07 


- 15.86 


—62 ± 16 


152 


—209 ± 46 


14.93 — 


15.49 


> 14.05 




—396 ± 115 


—919 ± 117 


22029224 


1.09 ± 0.19 


0.22 


14.84 


- 15.77 


—69 ± 16 


117 


— 146 ± 53 


14.81 — 


15.37 


14.37 - 


14.94 




—450 ± 39 


22036854 


2.08 ± 0.24 


0.24 


15.12 


- 15.81 


—91 ± 14 


58 


— 121 ± 37 


14.91 — 


15.81 


14.55 — 


15.30 


AAA 1 A f 

—444 ± 46 




22036912 


0.84 ± 0.21 


0.34 


14.73 


- 15.96 


—86 ± 16 


167 


—238 ± 87 


14.80 — 


15.29 


14.14 - 


14.49 


—297 ± 53 


—653 ± 58 


32010773 


0.63 ± 0.17 


0.19 


14.60 


- 15.70 


-197 ± 46 


45 


— 189 ± 126 


14.45 — 


15.03 


14.13 — 


14.60 


-511 ± 95 


—703 ± 43 


32011098 


0.59 ± 0.31 


0.33 


14.58 


- 15.94 


—89 ± 29 


87 


—78 ± 101 


14.54 — 


15.04 


14.56 — 


15.57 






32011099 


1.79 ± 0.25 


0.23 


15.06 


- 15.79 


-46 ± 13 


72 


-97 ± 24 


15.06 — 


16.12 


14.65 — 


15.30 


-356 ± 69 


-463 ± 105 


32011192 


1.85 ± 0.19 


0.22 


15.07 


- 15.77 


-154 ± 9 


97 


-154 ± 38 


14.96 - 


15.65 


14.74 _ 


15.30 


—527 ± 86 


—739 ± 74 


32011682 


1.30 ± 0.20 


0.23 


14.92 


- 15.79 


—85 ± 11 


96 


—246 ± 60 


14.90 - 


15.55 


14.18 - 


14.68 


—317 ± 71 


—521 ± 62 


32016857 


0.88 ± 0.18 


0.20 


14.75 


- 15.73 


-91 ± 15 


97 


—227 ± 82 


14.75 - 


15.32 


14.42 - 


15.53 


—297 ± 84 




32017112 


1.25 ± 0.25 


0.30 


14.90 


- 15.90 


—80 ± 18 


99 


—255 ± 106 


14.80 - 


15.38 


14.24 - 


14.88 






32022156 


0.95 ± 0.34 


0.43 


14.78 


- 16.06 


— 180 ± 48 


74 


—318 ± 182 


14.82 - 


15.49 


14.41 - 


14.98 


—355 ± 42 


—436 ± 27 


42006781 


2.07 ± 0.17 


0.15 


15.12 


- 15.60 


-70 ± 13 


59 


-143 ±28 


15.00 - 


16.11 


14.06 - 


14.42 


-565 ± 77 


-495 ± 68 


42006915 


1.29 ±0.10 


0.097 


14.92 


- 15.41 


-76 ± 10 


206 


-322 ± 75 


14.67- 


15.08 


13.91 - 


14.03 


-621 ± 117 


-576 ± 91 


42014101 


1.08 ± 0.35 


0.32 


14.84 


- 15.93 


-150 ±38 


107 


-281 ±251 


14.65 - 


15.17 


14.48 - 


15.29 


-514 ±68 


-480 ± 66 


42014154 


0.84 ±0.22 


0.25 


14.73 


- 15.82 


-105 ±30 


63 


-315 ± 163 


14.77- 


15.46 


13.93 - 


14.42 






42014585 


1.63 ±0.27 


0.27 


15.02 


- 15.86 


-177 ±30 


61 


-297 ± 85 


14.96 - 


15.94 


14.55 - 


15.20 


-579 ± 155 


-456 ± 47 


42014618 


1.29 ±0.14 


0.15 


14.92 


- 15.60 


-76 ± 13 


121 


-115 ±33 


14.83 - 


15.39 


14.43 - 


15.67 


-548 ± 107 




42014718 


2.06 ±0.31 


0.30 


15.12 


- 15.90 


-156 ± 19 


199 


-219 ± 77 


14.96 - 


15.49 


14.84 - 


15.90 


-712 ±216 


-651 ± 128 


42014732 


0.91 ±0.17 


0.16 


14.76 


- 15.63 


-81 ± 26 


79 


-52 ± 43 


14.59 - 


15.13 


14.47 - 


15.33 


-450 ± 82 




42021266 


1.06 ±0.20 


0.19 


14.83 


- 15.70 


-108 ±24 


76 


-388 ± 152 


14.75 - 


15.37 


14.28 - 


15.26 


-590 ± 163 




42022173 


0.81 ± 0.21 


0.23 


14.71 


- 15.79 


-134 ±27 


106 


-125 ±85 


14.67- 


15.20 


14.70 - 


15.85 






42022307 


1.15 ±0.09 


0.10 


14.87 


- 15.43 


-146 ± 10 


165 


-381 ± 79 


14.71 - 


15.18 


14.20 - 


14.66 


-554 ± 62 


-529 ±50 


42026243 




0.70 




16.27 


-302 ± 38 


93 


-476 ± 201 


15.17- 


16.22 


15.03 - 


16.29 




-843 ± 135 


42033991 


1.39 ±0.46 


0.40 


14.95 


- 16.03 


-139 ±41 


132 


-369 ±319 


14.69 - 


15.19 


14.52 - 


15.36 


-466 ± 136 





Note. — (1) DEEP2 survey identification. (2) Rest-frame equivalent width of Fe II 2374 absorption from direct integration of the trough. Error estimated from bootstrap resampling. 
(3) Upper limit (lcr) on Fc II 2261 equivalent width. Computed from an assumed lincwidth that is twice the FWHM of the line response function and the continuum S/N ratio. (4) 
Limits on Fe II column density for a fiducial covering factor of Cf — 1. For optically thin absorption, the linear portion of the curve-of-growth gives N(Fe II) = 6.40 x 10 14 cm W23745 
which we adopt as a lower bound. We derive an upper bound, JV(Fe II) = 9.06 x 10 15 cm 2 VK226ii by taking the 3cr upper limit on the Fe II 2261 equivalent width. (5) Doppler shift 
of the single velocity component model fitted to Fc II absorption as described in Section 3.3. (6) Doppler parameter of the systemic component in the two-component fits to the Fe II 
series. We adopt characteristic velocities of gas in the galaxy as probed by cooling radiation from H II regions. Specifically, we require the FWHM of the absorption optical depth profile 
match that of [OH] measured previously by Weiner etal. (2006), which gives b sys — \^2cr([OII]). (7) Doppler shift of the Doppler component fitted to Fe II lines in two-component 
model described in Section 3.3.2. (8) Limits on Fe II column density of ISM component in two-component fits. The model line profile for Fe II 2261 and 2374 were integrated and 
the above formulae, from comment #4, applied. (9) Limits on Fc II column density of Doppler component in two-component fits. The mode 1 line profile for Fe II 2261 and 2374 were 
integrated and the above formulae, from comment #4, applied. (10) Velocity of maximum Fc II 2374 absorption for J + 6(1) = 1. (11) Velocity of maximum Mg II 2796 absorption for 
I + S(I) = 1. 



TABLE 4 
Inflow Properties 



Object 


W(2374) 


W(2260) 


\og[N(Fe+)C f ] 


Vi 


biSM 


Vdop 


log[7V 


iSM(Fe+)C f ] 


log[N DO p(Fe+)C f ] 


Vmax (A2374) 


Vma*(A2382) 




(A) 


(A) 


( cm- 2 ) 


( km s- 1 ) 


( km s -1 ) 


( km s -1 ) 




( cm -2 ) 


( cm- 2 ) 


(km/s) 


(km/s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


12011428 


2.14 ±0.48 


0.44 


15.1 - 16.0 


154 ± 32 


86 


392 ± 124 


15.14 


-16.19 


14.71 -15.49 


416 ± 0.9 


523 ± 0.1 


12016019 


2.44 ±0.20 


0.19 


15.1 - 15.7 


48 ± 14 


97 


258 ± 40 


15.17 


-16.15 


> 13.58 


655 ± 34 


723 ± 0.6 


12019962 


1.88 ±0.34 


0.32 


15.0 - 15.9 


104 ± 31 


113 


122 ±97 


14.68 


-15.20 


13.86 -15.15 


621 ± 9 


472 ± 1 


12019996 


1.69 ±0.31 


0.30 


15.0 - 15.9 


117 ± 24 


158 


386 ± 138 


14.92 


-15.47 


14.56 -15.57 


510 ± 3 


577 ±5 


12100420 


1.86 ±0.28 


0.26 


15.0 - 15.8 


71 ± 18 


109 


190 ± 45 


15.14 


-15.98 


14.37 -14.82 


539 ± 1 


641 ± 1 


22004858 


1.94 ±0.14 


0.14 


15.0 - 15.6 


61 ± 11 


151 


71 ± 18 


15.03 


-15.63 


14.76 -15.53 




737 ±2 


22005270 


0.98 ±0.23 


0.25 


14.8 - 15.8 


401 ± 43 


145 


505 ± 588 


14.39 


-14.71 


14.37 -14.81 


546 ± 2 


724 ±9 


32020468 


2.36 ±0.69 


0.62 


15.1 - 16.2 


123 ± 39 


95 


164 ± 79 


15.21 


-16.33 


14.82 -15.97 


465 ± 1 


311 ± 3 


42006875 


1.01 ±0.12 


0.12 


14.8 - 15.5 


51 ± 14 


62 


256 ± 50 


14.87 


-15.66 


13.86 -14.21 




663 ±3 



Note. — (1) DEEP2 survey identification. (2) Rcst-framc equivalent width of Fc II 2374 absorption from direct integration of the trough. Error estimated from bootstrap resampling. 
(3) Upper limit (lcr) on Fe II 2261 equivalent width. Computed from an assumed lincwidth that is twice the FWHM of the line response function and the continuum S/N ratio. (4) 
Limits on Fe II column density for a fiducial covering factor of C / — 1. For optically thin absorption, the linear portion of the curve-of-growth gives 7V(Fe II) — 6.40 X 10 14 cm - 2 VK2374, 
which we adopt as a lower bound. We derive an upper bound, Af(Fe II) = 9.06 X 10 15 cm from the 3cr upper limit on the Fe II 2261 equivalent width. (5) Dopplcr shift of the 

single velocity component model fitted to Fc II absorption as described in Section 3.3. (6) Doppler parameter of the systemic component in the two-component fits to the Fe II series. 
We adopt characteristic velocities of gas in the galaxy as probed by cooling radiation from H II regions. Specifically, we require the FWHM of the absorption optical depth profile 
match that of [Oil] measured previously by Wcincr ct al. (2006), which gives b sys — y/2a([OII]). (7) Dopplcr shift of the Doppler component fitted to Fe II lines in two-component 
model described in Section 3.3.2. (8) Limits on Fe II column density of ISM component in two-component fits. The mode 1 line profile for Fe II 2261 and 2374 were integrated and 
the above formulae, from comment #4, applied. (9) Limits on Fe II column density of Doppler component in two-component fits. The model line profile for Fe II 2261 and 2374 were 
integrated and the above formulae, from comment #4, applied. (10) Velocity of maximum Fe II 2374 absorption for J + 5(1) = 1. (11) Velocity of maximum Mg II 2796 absorption 
for / + <5(J) = 1. 
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TABLE 5 
Properties of Composite Spectra 



Feature Vi (Fe II) EW(A2374) V m ax(MglT) EW(A2796) 



Spectrum 


( km s 1 ) 




(A) 


( km s 


- 1 ) 




(A) 




(1) 


(2) 




(3) 


(4) 






(5) 




ijOwcst iviass i^o.oo ^ logi^-M* / Mq j <. y.cso j 




1 


.ZD 


fi/IQ J. 

— 04y it 


OD 


A 

u. 


01 ± u. 


ni 
Ul 


Intermediate Mass (y.oo S logi^iVZ*/ Mq j < lU.ooj 


—32 i 6 


1 


.47 


— 770 ± 


103 


1, 


.40 ± 


03 


xligncst Mass (^lU.oo ^ log^iVi*/ Mq J ll.zyj 


— ZO ± 


1 

1 




— / lo ± 


iin 
lly 


1 

1. 


.Do ± U 


-UO 


Bluest (0.189 < U — B < 0.519) 


-32 ±8 


1 


.31 


-693 ± 


86 





.72 ± 0. 


.02 


Bluer (0.519 < U - B < 0.662) 


-27 ± 6 


1 


.48 


-649 ± 


91 


i 


.17 ±0. 


02 


Blue (0.662 < U — B < 0.907) 


-26 ± 7 


1 


.52 


-888 ± 


173 


i 


.85 ± 0. 


.07 


Lowest Luminosity (-17.67 > M s — 5 log h > -20.03) 


-23 ±9 


1 


.19 


-666 ± 


88 





72 ±0. 


02 


Intermediate Luminosity (-20.03 > M B - 5 log ft > -20.79) 


-18 ±7 


1 


.56 


-667 ± 


88 


1 


.08 ± 0. 


.01 


Highest Luminosity (-20.79 > M B - 5 log ft > -22.19) 


-39 ±6 


1 


.65 


-782 ± 


92 


1 


.79 ±0. 


.03 


Lowest SFR (0.86 < 5FR( M yr" 1 ) < 12.24) 


26 ± 10 


1 


.48 


-629 ± 


104 





.81 ±0 


.04 


Highest SFR (12.24 < SFR{ M© yr" 1 ) < 97.81) 


-38 ± 8 


1 


.47 


-755 ± 


91 


1 


.88 ±0 


.02 


Lowest sSFR (1.00 < sSFR~ 1 (Gyr) < 15.08) 


19 ± 10 


1 


.55 


-779 ± 


136 


1 


.59 ±0 


.06 


Highest sSFR (0.178 < sSFR^iGyr) < 1.00) 


-24 ± 9 


1 


.41 


-610 ± 


81 





.92 ±0 


.02 


3ct Fell Outflow 


-119±6 


1 


.47 


-901 ± 


99 


1 


.59 ±0 


.02 


3cr Fell Inflow 


95 ±8 


1 


.89 


-559 ± 


113 


1 


.25 ±0 


.07 



Note. — Col. 1 - Subsamplc averaged together to form composite spectrum. Col. 2 — Dopplcr shift of the Fc II absorption 
measured as described in Section 3.3. Col. 3 — Equivalent width of the Fe II A2374 absorption trough. Col. 4 — Maximum 
blucshift, Vmax, of Mg II A2796 absorption measured as described in Section 3.4. Col. 5 - Equivalent width of the Mg II A2796 
trough from to Vmax- 
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TABLE 6 

Correlations Among Properties of Individual Galaxies 



Vo£ UeLIllrl by 




vv A2374 


v * \2796\ v ^ U J 


Vmax 


Mgll 


M, 


0.099 (1.21) a 


0.139 (1.70) a 


0.320 (1.87) c 


-0.305 (1.78) c 


0.493 (2.9) c 


M B - 5 log h 


0.048 (0.58) a 


-0.133 (1.62) a 


-0.338 (1.97) c 


0.097 (0.56) c 


-0.495 (2.9) c 


U - B 


0.123 (1.50) a 


0.116 (1.42) a 


0.443 (2.58) c 


-0.205 (1.20) c 


0.685 (4.0) c 


z 


-0.054 (0.66) a 
0.021 (0.26) a 


0.180 (2.20) a 


0.066 (0.38) c 


0.143 (0.84) c 


0.401 (2.34) c 


SNR 


-0.065 (0.80) a 


-0.165 (0.96) c 


-0.059 (0.34) c 


-0.050 (0.29) c 


Vi(Fell) 




0.221 (2.70) a 


-0.1421 (0.83) c 
0.451 (2.63) c 


0.160 (0.93) c 
-0.033 (0.19) c 


0.159 (0.92) c 
0.410 (2.39) c 




0.221 (2.70) a 




V m ax(Mg II) 


0.160 (0.93) c 


-0.033 (0.19) c 


-0.451 (2.63) c 




-0.302 (1.76) c 


W X2796 (V < 0) 


-0.142 (0.83) c 


0.451 (2.63) c 




-0.451 (2.63) c 


0.717 (4.2) c 


SFR 


-0.225 (1.40) b 


-0.005 (0.03) b 


0.086 (0.19) d 


0.257 (0.57) d 


0.714 (1.60) d 


sSFR 


-0.257 (1.61) b 


-0.079 (0.49) b 


-0.257 (0.58) d 


0.086 (0.19) d 


0.314 (0.70) d 



NOTE. — Wc computed the Spearman rank-order correlation coefficient, r$, between the continuum S/N 
ratio near Fc II and various galactic properties. To determine whether a significant correlation was present, 
we examined the variance in the sum-squared difference of ranks, D, which is approximately normally 
distributed as described in chapter 14 of Press et al. (1992). The value in parentheses is the number of 
standard deviations from the null hypothesis (no correlation between the two quantities). 



a Samplc of 150 blue-cloud galaxies with measured Fell Dopplcr shift Vi . 

b Sample of 40 blue-clouds galaxies with measured Fell Doppler shift Vi and extinction-corrected, ultravi- 
olet SFR. 

c Samplc of 35 high S/N ratio spectra with measured maximum Mg II blucshift y max . 

d Sample of 6 high S/N ratio spectra with measured V^a^fMgll) and extinction-corrected, ultraviolet 
SFR. 



